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Progress in Dynamic Stability and Control 
Research 


WILLIAM F. MILLIKEN, JR.* 


Cornell Aeronautical Laboratory 


SUMMARY 


This report presents a discussion of an extensive flight and 
analytical research program in dynamic stability and control 
which is being performed at the Cornell Aeronautical Laboratory. 
The basic objective of this program is that of bringing dynamic 
analysis closer to the needs of the airplane, missile, and auto- 
matic control designer, through a simplification of approach and 
by the development of methods for obtaining quantitative dy- 
namic information in flight. Automatic control, in the form of 
modified automatic pilots, is used throughout the flight program. 

In Part I of this report a review of the significance, background, 
and scope of the work is given. The important technical as- 
pects are then discussed in Parts II and III, which survey the 
work done in the longitudinal motion over the last 2'/: years, 
while in Part IV the future direction of the overall program is 
outlined. Finally, an opportunity is taken to summarize in the 
Appendix other dynamic research, both as a background for the 
present efforts and because of its inherent interest at the present 


‘time. 


The program to date has resulted in the development of a num- 
ber of flight and analytical methods for the measurement and 
handling of dynamic phenomena. Some of these are immediately 
applicable to design problems. The most outstanding accomplish- 
ment of the program, from an overall point of view, has been its 
influence in furthering the use of control forcing functions (and 
automatic control as a means of applying them) in the determi- 
nation of aerodynamic information. 


INTRODUCTION 


ee ARE CERTAIN PROBLEMS in the design of air- 
planes and other types of flying devices which de- 
pend for solution upon a detailed knowledge of the 
transient velocity conditions that result from flight 
control use and atmospheric or other disturbance. As- 
sociated with the time history of the motion—unsus- 
ceptible to rigorous solutions by static analysis alone— 
these problems come under the broad classification of 
dynamic stability and control. 


Presented at the Flight Test Session, Fifteenth Annual Meet- 
ing, I.A.S., New York, January 28-30, 1947. Revised and re- 
ceived May, 1947. 

* Flight Research Department. 


Rapid developments in missiles and conventional air- 
craft have placed a premium on new methods and tech- 
niques that will enable the potential utility of dynamic 
analysis to be realized in design. 

The present report reviews the salient features of a 
full-scale dynamic stability and control research pro- 
gram now in progress at the Cornell Aeronautical Lab- 
oratory and its relation to allied investigations. Dur- 
ing this program extensive measurements are being 
made of airplane response to control motions that are 
applied by means of automatic control. 

In the first part of this report an examination is made 
of the general background and objectives of the work, 
while the remaining parts are devoted to specific 
phases of the program. 


Part I—GENERAL 


The Necessity for Dynamic ‘Research 


While, from the earliest days of the airplane, it was 
appreciated that a complete understanding of stability 
and control involved dynamic analysis (taking into ac- 
count not only the static attitude restoring terms but 
all of the aerodynamic and inertia effects arising from 
the motion), relatively little was accomplished in the 
ensuing years to make this a commonly accepted tool 
in design. In general, little was done toward simplifica- 
tion of the mathematical approach, determination of 
the operating limits within which the classical assump- 
tions were valid, or devising methods of obtaining 
specific experimental data (either from wind-tunnel or 
flight test) for use in design calculations and full-scale 
changes and improvements. This lack of engineering 
adoption is, in part, explainable by the fact that it was 
possible to develop satisfactory airplanes on the basis 
of evolutionary design experience, static wind-tunnel 
test, and cut-and-try modification stemming from flight 
tests on the prototype and service use on subsequent 


493 


e 
‘ 
; 
DF 
4 
4 
| 
JER 


494 JOURNAL OF THE AERONAUTICAL SCIENCES—SEPTEMBER, 


articles, while, in addition, the existing specifications 
at that stage of the airplane art did not, in most cases, 
warrant the cost of seeking more exact design data. 
However, vast increases in airplane, missile, and pilotless 
aircraft performance and new uses of automatic control 
are making successful and economic design possible only 
through a detailed knowledge of dynamic behavior. 

An examination of current missile development and 
that slated for the forseeable future serves to emphasize 
the wide range of specialized problems that require both 
a knowledge of aerodynamic response characteristics 
and a know-how in their application, which have hardly 
been a part of the traditional aerodynamicist’s back- 
ground. For example, the ability to design guidance 
equipment for beam- and radar-controlled missiles pre- 
supposes sufficient information to perform a closed loop 
analysis in which the component response characteris- 
tics are synthesized for the overall system behavior. 
This information, while readily available for the elec- 
tronic elements of a controlling system, has yet to be 
provided satisfactorily for the aerodynamic ones. In 
these applications the aerodynamicist is being forced to 
face such fundamental issues as the time constants re- 
lating flight path and attitude changes to control dis- 
placement and the result of application of various con- 
trol forcing functions. He is being required to gear-up, 
from his essentially static aerodynamics, to such con- 
cepts as feedback and complex transfer functions that 
are used in dynamic electrical and servomechanism de- 
sign. 

In the man-carrying aircraft field the use of dynamic 
analysis is equally urgent. As automatic control daily 
finds new uses, its designers stress the need for more 
aerodynamic response information. Improvement in 
the rough-air riding qualities and safety of transport 
aircraft when flown on automatic control will certainly 
be demanded in the near future, along with advances in 
push-button automatic flight control methods. Radar 
bombing and tracking devices must be considered in- 
herently with the dynamics and automatic control of 
the aircraft in which they are installed. A corollary to 
the design of aircraft whose performance approaches or 
extends into the trans- and supersonic régimes is that 
many features of the external configuration (such as 
sweepback, etc.) are fixed by the performance require- 
ments themselves, thus enhancing the difficulties of the 
stability and control designer and necessitating more 
precise methods of analysis. Formidable dynamic 
problems obviously face the designer in the transonic 
region and at the stall. Finally, in the case of the large 
or unconventional aircraft, the prototype investment 
may become so huge that, in the light of increasingly 
stringent flying qualities specifications, it will be no 
longer economical (or acceptable to the operating 
agencies) to by-pass a complete dynamic design study— 
a trend that is already reflected in the incorporation of 
dynamic requirements in military and civilian flying 
characteristics specifications. 
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In short, the manufacturers of either military flying de- 
vices or commercial airéraft will, to an increasing degree, 
be called upon to predict and demonstrate quantitatively 


. the dynamic behavior of their products. On the one hand, 


it is an inescapable fact that the aircraft is a dynamic 
machine whose more precise examination is an outgrowth 
of its widening utility; on the other, it is equally evident 
that the aerodynamic designer is still not equipped to solve 
dynamic problems with the desired degree of facility. It is 
in relation to these facts that the necessity and opportunity 
for new approaches in dynamic research lie. 


Initiation of a Research Program 


Two and one-half years ago, from observation of these 
trends, a project was initiated by the author at the 
Cornell Aeronautical Laboratory (then the Curtiss- 
Wright Research Laboratory) as an entrée to a broad 
research in the field of dynamic stability and control and 
its application to future flying devices. The overall re- 
search, comprising several projects at the present time, 
and a major endeavor of the Laboratory, has the general 
objective of furthering the use of dynamic analysis by 
the aerodynamic designer through a simplification of 
approach and by the provision of full-scale numerical 
data. Because of the limited familiarity that most 
aeronautical engineers have with dynamic stability and 
control, the project seemed particularly suited to the 
long-range objective of the Laboratory, which is one of 
combining research and education. Up to the present 
time, the research has been entirely performed in the 
Laboratory’s Flight Research Department, which has 
been, in many respects, organizationally tailored to the 
unique requirements of dynamic research in flight. 

The program has enjoyed, in addition to the early 
support of Curtiss-Wright, the primary sponsorship of 


the Aircraft Laboratory, Engineering Division, Air: 


Materiel Command, which has contracted for a major 
block of the experimental and theoretical work now 
under way. Currently, a portion of the research (in- 
cluding an application of the flight oscillation technique 
to a particular design problem) is being performed for 
the Instrumentation Laboratory at Massachusetts In- 
stitute of Technology, while the early exploratory phase 
of the program was done on a mutually cooperative 
basis with the Flight Research Branch, Flight Test 
Division, Air Materiel Command, and the Minneapolis- 
Honeywell Regulator Company. The Sperry Gyroscope 
Company, Inc., has extended informal cooperation 
during a current phase of the program which utilizes 
their A-12 type of automatic pilot. 

It is the purpose of the present report to deal with 
the general progress of this research, viewed as a whole 
rather than as individual projects. Being in the nature 
of a survey, no attempt is made to present more than 
typical samples of the large amount of flight data that 
have been collected. Theoretical developments that 
have progressed in parallel with the flight program are 
summarized only with respect to the composite picture, 
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since it is hoped that detailed presentations will be 
made by their respective authors at an early date. 
This report is published with the full cognizance of, 
and appreciation for, the interest of the various spon- 
soring agencies. 


Ideas Underlying the Approach 


In initiating a dynamic research program it was, of 
course, appreciated that a great deal had already been 
accomplished in the field and that many different 
avenues of possibility were competently being explored. 
To avoid an overlapping of effort and ensure an orienta- 
tion of objectives, a survey was made of previous dy- 
namic research.* An outline of the proposed program! 
was then circulated to interested automatic pilot and 
aircraft manufacturers, the National Advisory Com- 
mittee for Aeronautics, and others for their suggestions. 

In brief, it was apparent that aircraft dynamic sta- 
bility and control considerations had for years (and 
almost without exception) revolved about the so-called 
classical equations of motion.? In these equations, 
which rigorously hold only for small deviations from 
steady symmetrical flight, the airplane was taken as a 
rigid body acted upon by aerodynamic, inertia, and 
gravitational forces and moments. The unbalanced 
aerodynamic forces and moments were expressed as 
the sum of several elements, each of which varied lin- 
early with an increment in a respective linear or angular 
velocity component, these force and moment curve 
slopes being referred to as the “‘stability derivatives.”’ 
The total forces and moments were largely assumed 
independent of acceleration components. 

In considering further research it was immediately 
apparent that additional evidence, preferably of a full- 
scale experimental nature, was needed to substantiate 
the assumed variation of the aerodynamic forces and 
moments on modern aircraft. It was later realized that 
the constancy or the unique existence of the simple 
velocity derivative was questionable in the more rapid 
modes of motion, which meant that the actual condi- 
tions under which the various types of derivatives 
would adequately describe the motion were largely un- 
explored. Furthermore, few measurements had been 
made in flight to substantiate even the simple individual 
velocity derivatives as determined from the wind tun- 
nel, largely for the reason that no generally accepted 
flight-test method had been developed (although a 
considerable amount of closely allied data, principally 
static, were obtainable from routine flying qualities 
tests). Such substantiation, it was felt, might go far 
toward suggesting new wind-tunnel techniques or the 
clarification of design criteria for a new type of dy- 
namic wind-tunnel balance; while, on the other hand, 
a simple scheme for measuring dynamic characteristics 


* Because of its importance to an understanding of the formula- 
tion of objectives for the present program, an up-to-date synop- 
sis of this survey is given in the Appendix. 


in flight would be of outstanding significance for full- 
scale research at Mach Numbers well beyond the range 
of dynamic or reliable static wind-tunnel data (it being 
evident that a complete dynamic determination would 
include all of the restoring derivatives that correspond 
to the static stabilities in common use). 

One was also led to inquire, in the light of experience 
in other branches of science (notably electrical, where 
semimathematical or graphical methods are in use for 


' handling dynamic phenomena), whether there was not 


some simpler means of expressing the dynamics of the 
aircraft—at least for certain problems—which would 
stress the physical aspects, avoid the cumbersome 
mathematics attendant on the use of the equations with 
control moments present, and, in general, minimize the 
limitations of the classical approach. It seemed de- 
sirable to tap the broad experience of the electrical and 
electronics engineers, who are the recognized experts 
in handling dynamic phenomena. Although most 
aeronautical designers felt in a vague way that an air- 
craft was a ‘‘spring’’ system, they had little familiarity 
with its characteristics even for the more important 
practical modes. 

Finally, it was concluded that almost complete aerody- 
namic information was lacking which would allow the 
engineer to match automatic control and aircraft charac- 
teristics, a subject that could well be the most important 
phase of any dynamic investigation. 

The above considerations led to the decision that the 
program should be organized as a full-scale research in 
flight, backed by closely coordinated analytical inves- 
tigations. Among the advantages of performing the 
research in flight was that of being able to observe and 
experience physically the phenomena under study. (In 
the author’s opinion, ‘‘on-the-spot’’ engineering intelli- 
gence is a flight research resource too often overlooked 
and one, of incalculable value in the dynamics field 
where a consolidation of theory and practice is so im- 
portant). Jt was believed that a thorough understanding 
of the dynamics of the airplane—which because of its size, 
flexibility, and speed, is so relatively easy to study and 
make multiple measurements on—is a first requisite to 
unraveling the dynamic behavior of the missile. 

It was also decided to use automatic control; first, 
in the capacity of a flight-test tool (an application orig- 
inally suggested by Kutzler and Johnston*), which 
would allow the airplane to be flown repeatedly through 
precise specialized maneuvers; and, second, for the 
excellent opportunity presented for quantitatively 
studying its behavior in relation to the aerodynamics of 
the airplane. 

With automatic control present, the most promising 
method of test appeared to be that of applying known 
(predetermined) control motions—i.e., control forcing 
functions—simultaneously recording the airplane re- 
sponse. For some reason this well-known and powerful 
method of handling dynamic phenomena in general 
had seen but little aeronautical application. 
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In summary then, it was decided to do the entire job in 
flight, taking the fullest advantage of whatever automatic 
control could offer in attacking the problem from the forcing 
function standpoint. While it was hoped that utilizing 
electrical engineering experience would suggest ways and 
means of improving upon the classical theory, attention 
would first be directed to the determination of the velocity 
derivatives and the establishment of valid operating 
limits. 


Viewed as an integrated whole, the above approach ' 


differed markedly from that of other dynamic research; 
the additional insurance, it was hoped, that the work 
would be complementary to allied endeavor. 


Summary of Specific Objectives 


On the basis of the ideas outlined above, it was pos- 
sible to state a set of definite long-range objectives. 
These objectives are as follows: 

(1) To develop automatic control and instrumen- 
tation methods and techniques for measuring air- 
plane dynamic stability and control characteristics in 
flight. 

(2) To perform comparative tests with existing 
types of automatic pilots to ascertain characteristics 
most suitable for their use in item (1); to suggest and 
test modifications enhancing their use in that respect. 

(3) To amass a fund of numerical forcing function 
response data from flight over a wide range of control 
application rates and for various aircraft and operating 
conditions. 

(4) To study the relationship between flight-test 
response information and the classical derivatives; to 
compare the response information, properly reduced, to 
wind-tunnel and theoretical information to further 
confirm or establish limits of its applicability. 

(5) To suggest new and simpler theoretical ap- 
proaches, or wind-tunnel techniques, falling in line 
with the flight experience and easily usable in practical 
design. 

(6) To study and demonstrate the feasibility of 
using the experimental response methods in high Mach 
Number man-carrying airplane applications and in 
missile and pilotless aircraft design. 

(7) To enhance the measurement and application 
of airplane dynamic stability and control information 
for use in the design of automatic control, servomech- 
anisms, and power boost. 

(8) To provide dynamic information in a form 
suitable for simulator design or for checking the per- 
formance of simulators. 

(9) To adapt the dynamic approaches used in 
other branches of science to the needs and problems of 
the aircraft designer. 

(10) Tostudy and research experimentally upon the 


application of automatic control for unstable and un- 
conventional aircraft, stall control, ride control, and 
other specialized flight problems. 
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(11) To study experimentally, if such appears prac- 
tical, the combined longitudinal and lateral motion. 

(12) To simplify such new automatic control meth- 
ods as are developed to the point where they may be 
incorporated into routine flight testing (as opposed to 
their pure research use); to devise and perfect flight 
methods for measuring dynamic information which 
may be applied independently of automatic control. 

Many secondary objectives outside the scope of this 
report were also established. 


Organization and Remarks 


During the last 2'/: years, the effort expended on this 
research has grown to about 60 per cent of that of the 
Flight Research Department, an organization of 70 
people. In addition to the flight crew, man power is 
available for project, instrumentation, mechanical, 
photographic, and theoretical and analysis engineering 
functions. The Flight Research Department is fortu- 
nate in possessing all of those organizational elements 
necessary to a research program in flight and an as- 
sociated ground analytical investigation in a closely in- 
tegrated, singly directed organization. Of even more 
consequence to the present research is the fact that in 
one relatively small group resides specialized experience 
in aerodynamics, electronics, and applied mathematics, 
as well as in flight testing. More than 250 hours of 
flying have been devoted to the dynamic stability and 
control work since its inception, the greatest portion of 
this flying being on service airplanes of the medium 
bomber category. 

Among the problems faced in a research of this na-, 
ture are those associated with a stacking of the well- 
known (but often minimized) intricacies of multiple- 
channel dynamic recording equipment with those of 
adapting automatic control to new uses. To these there 
must be added the usual and numerous routine prob- 
lems of carefully controlled flight testing. Obviously, 
great care must be taken in the measurements so that 
the dynamics of the airplane (if under test) are iso- 
lated from the dynamics of the automatic control and/ 
or the recording equipment. Brief discussions of 
some of the interesting instrumentation problems 
encountered are included in later sections of this re- 
port. 

The research as a whole has attained many of the 
goals initially set and has thoroughly established the 
longitudinal oscillation forcing function technique for 
obtaining dynamic stability and control data. It is 
now progressing upon an investigation of the step func- 
tion approach, the application of the longitudinal tech- 
nique to particular problems, and the oscillation tech- 
nique to the lateral and combined lateral and longitu- 
dinal motions. While the future direction of the work is 
discussed in a later section of this report, it may be 
mentioned here that the program bids well to be a con- 
tinuing endeavor of increasing magnitude. 


P 


Intro 


In 
descr 
(abot 
nami 
plane 
Four’ 
pared 
the 
prese 

Th 
25J 
C-1t 
and 
spons 
mal a 
move 
be in 
merel 
forcin 
by a. 
the n 
prove 
contre 
impos 
thus « 
path. 
affect 
meast 
reasol 

The 
forcin 
least 
thoug 
time 
deter1 
contr« 
functi 
on the 
at lea 
worki 
on dis 
low n 
reasor 

The 
need 
chani 


menta 
sign. 
of thi 
desig 


anis 
Instrq 


Fo 
sirab 


| this 
the 
f 70 
er is 
‘ical, 
ring 
rtu- 
ents 
| as- 
y in- 
nore 
it in 
ence 
tics, 
s of 
and 
n of 
ium 


DYNAMIC STABILITY AND CONTROL RESEARCH 497 


Part II—LONGITUDINAL OSCILLATION METHOD 


Introduction 


In accordance with the general method of approach 
described in Part I, a series of tests were initiated 
(about 2 years ago) for the determination of the dy- 
namic longitudinal response characteristics of an air- 
plane. The results of these tests were reported‘ at the 
Fourteenth Annual I.A.S. Meeting but were never pre- 
pared for publication. Acknowledgment is made to 
the coauthors for the use of a part of this material in the 
present report. 

The tests were performed on a North American B- 
25] airplane equipped with a Minneapolis-Honeywell 
C-1 type automatic pilot, and the longitudinal stability 
and control characteristics were obtained from the re- 
sponse of the airplane, as measured by variation in nor- 
mal acceleration and pitch angle when the elevator was 
moved with a continuous sinusoidal motion. It will 
be immediately recognized that the test method is 
merely one of applying a sinusoidal elevator control 
forcing function and determining the airplane response 
by a measurement of the most convenient variables of 
the motion. The automatic pilot, suitably modified, 
proved to be a logical device'for applying this specialized 
control motion, which, it should be emphasized, is super- 
imposed upon the basic function of the automatic pilot, 
thus enabling the maintenance of a desired mean flight 
path. This superposition has been shown to in no wise 
affect the results of the tests (which are based on 
measured elevator angle) if the air conditions are 
reasonably smooth. 

The application of a sinusoidal elevator control 
forcing function had previously been explored by at 
least one automatic pilot manufacturer (Sperry), 
though the present tests probably represent the first 
time that the technique has been used in the specific 
determination of full-scale aerodynamic stability and 
control data. The decision to use the sinusoidal forcing 
function, in preference to some other type, was based 
on the desirability of obtaining (during the initial tests, 
at least) steady-state information. The advantages of 
working with the longitudinal motion were emphasized 
on discussion with Koppen, who predicted the relatively 
low normal acceleration that was obtained for any 
reasonable elevator motion. 

The potentialities of the sinusoidal forcing function 
need hardly be mentioned. Through its use in me- 
chanical and electrical systems, a large fund of experi- 
mental information has been collected for rational de- 
sign. While there was some doubt as to the reduction 
of this type of information for direct use in airplane 
design, its utility in automatic control and servomech- 
anism adaptation has never been questioned. 


Instrumentation, Procedure, and Tests 


For the purpose of developing the method it was de- 
sirable to obtain a continuous time history of all of the 


variables, while for routine use more preferable means 
of obtaining the pertinent data may be evolved. For 
the present tests continuous recordings of the following 
items were obtained on an oscillograph: (1) elevator 
angle; (2) elevator hinge moment; (3) pitch angle; 
(4) normal acceleration at nose, center of gravity, and 
tail; (5) angle of attack; [(6) bending stresses in hori- 
zontal tail spars]. 

The large gravitationally controlled flight gyro of the 
automatic pilot itself was used as the source of the 
pitch angle signal, while Consolidated reluctance-type 
accelerometers were used for the normal accelerations. 
The instrumentation was so installed in the compart- 
ment adjacent to, and immediately behind, the pilot's 
location that the operator could also control the auto- 
matic pilot and, without moving, observe the open- 
type photorecorder, the oscillographic equipment, the 
pilot’s instrument panel, and, in addition, see forward 
and to the side through the copilot’s windows. By this 
strategic location of ‘‘intelligence,’”’ a much clearer 
physical picture was obtained than from the mere ex- 
amination of recorded data. Fig. 1 shows a general 
view of the instrumentation and other equipment, 
while Fig. 2a is a typical record from the oscillograph. 
For these tests no filters for the reduction of hash were 
installed in the, oscillograph channels, although they 
were later introduced (the final development being 
shown in Fig. 2b). The elevator motion, obtained by 
the adaptation of an electronic oscillator control to the 
automatic pilot, as closely approximated the sinusoidal 
as could be expected of the basic constant speed servo- 
motor design. Fig. 3 shows a typical time history of 
the pitch angle, c.g. acceleration, elevator angle, and 
hinge moment variation. 

The test procedure consisted of stabilizing at a given 
flight condition and establishing an oscillation of small 
amplitude (approximately +11/;° elevator angle) at a 
given frequency. When steady-state conditions were 
attained, simultaneous oscillographic and photo-ob- 
server records were made for at least five complete os- 
cillation cycles. This procedure was repeated for the 
entire band of frequencies, ranging in period from 1 to 7 
sec. and giving approximately 25 points on the response 
curve for a given flight condition. Tests were performed 
in level flight, climb, and descent, but for the present 
purposes only the level flight runs at equivalent air 
speeds* of from 135 to 215 m.p.h. will be discussed. 
(Excellent data have, however, been obtained in climbs 
and descents.) 

It was found that the airplane oscillated steadily 
with no sensible change in air speed or longitudinal accel- 
eration and that for any reasonable value of elevator 
amplitude (< +3°) the normal acceleration increment 
at the center of gravity was less than +0.6g. Because 
of the combined rotation and vertical translation, the 
apparent center of rotation of the airplane, as seen in 


* Equivalent air speed is defined as the true air speed multi- 
plied by the square root of the density ratio. 
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AUTOPILOT 
CONTROL PANEL 


AUTOPILOT 
JUNCTION BOX} 


INSTRUMENTATION 
SWITCH PANEL 


General view of instrumentation, B-25J airplane. 


Fie. 1. 


Early records 


Typical oscillograph records from flight. 
with no electrical filter. 


Fic. 2a. 


reference to the horizon, was about 6 ft. ahead of the 
center of gravity and moved progressively forward with 
increases in air speed. While these figures are to some 
extent functions of airplane size, speed, etc., they serve 
to indicate that the method is physically practical. 
When the airplane was oscillated for a sufficiently 
long period, an approach to the phugoid mode was 
noticeable in that the air speed and altitude no longer 
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RECORDING 
OSCILLOGRAPH 


| 


Fic. 2b. Typical oscillograph records from flight. Recent 


records with electrical filters. 


remained exactly constant. (Further discussion of this 
point is presented later.) 

Experience with the method has indicated that the 
instrumentation and procedure are basically satisfac- 
tory and that, in general, it is a technique satisfactory 
throughout the speed range and in any attitude nor- 
mally experienced in this class of airplane. The test 
scatter of the pitch angle data was later reduced by 
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Fic. 3. Experimental time history of elevator angle, pitch 
angle, resultant normal acceleration of c.g., and elevator hinge 
moment for B-25J airplane at 175 m.p.h. equivalent air speed 
w; = 2.36 rad. per sec. 
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Fic. 4. Experimental variation of pitch angle amplitude with 
frequency of elevator oscillation for the B-25J airplane. 


suitable modification to the instrumentation, improving 
the readability of the data from the oscillograph record. 


Results 


Basic experimental results from these tests are shown 
in Figs. 4 and 5, where the magnitudes of the pitch (or 
attitude) angle (@, in degrees), and the resultant normal 
acceleration at the center of gravity (m, in g units), is 
in each case divided by the magnitude of the elevator 
angle (6, in degrees) and plotted vs. frequency in radians 
per second. These data are all reduced from time his- 
tories (as of Fig. 3) by halving measured double ampli- 
tudes averaged over a few cycles of the motion. The 
scatter of the test points, particularly in the case of the 
pitch angle, has since been reduced as experience in the 
instrumentation and data reduction has been ac- 
quired. Unfortunately, in these preliminary tests, the 
scatter did not permit accurate phase-angle measure- 
ment for the tests under discussion. 


Analysis and Discussion 


During the tests it had been observed that, through- 
out the frequency range covered, no noticeable change 
occurred in the air-speed indicator reading. Although 
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Fic. 5. Experimental variation of resultant normal accelera- 
tion of center of gravity with frequency of elevator oscillation for 
B-25J airplane. 


air-speed systems are notorious for lag, it was neverthe- 
less felt that the motion was occurring at substantially 
constant speed (an assumption that has since been 
substantiated is discussed further in Part III). It was 
then possible to make a simplified theoretical compari- 
son with the test data by eliminating the X equation 
and such terms as would arise from speed variation in 
the remaining classical longitudinal stability equations. 
If the dimensional equations* of reference 2 are thus 
modified, and the elevator forcing function, 5M; X 
sin wt, introduced (on the assumption that the forcing 
function is entirely one of pitching moment), a solution 
may be obtained by straightforward methods for the 
amplitudes of the pitch angle and acceleration response; 
thus, 


+ Gale? 
5 + (C2 — = 
Uw 


= M, 
+ (C2 — + 


where 


C. = + M, + 
= ZyM, — MyU 


— = |n|/s 


U = steady forward velocity 

n = resultant normal acceleration as 
measured by accelerometer at air- 
plane center of gravity 


6 = amplitude of elevator angle 
0 = angle of longitudinal airplane refer- 
ence axis to horizontal 
w = velocity increment along the Z axis 
M;, Zw, Mz = dimensional stability derivatives 


* For many comparative computations the use of the dimen- 
sional equations has been found convenient. 
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TASLE 


VALUES OF PARAMETER USED IN CALCULATION 
OF THEORETICAL RESPONSE CHARACTERISTICS 


E.AS.= 175 M.P.H. 


P = 0.00165 SLUGS/FT? 
U = 308 FT./SEC. 


= 0.0816 


(4 = 0.0769 


$= 610 FT? = -0,.0049 
9.67 FT. 

132 FT? (3%) = -0.017 

25.9 FT. (ss 0.45 

m= 840 SLUGS 7s 


Iy = 50,000 SLUG FT: hee 1.0 (POWER ON) 


(24), EUS -o864 i/sec. 
PUSE -9.00842  1/FT. SEC. 


t 


t Mq 4.57 FT./SEC 


2 


. sec? 


Numerical values of the various parameters calcu- 
lated from the best wind-tunnel and theoretical data 
available and for an equivalent air speed of 175 m.p.h. 
are summarized in Table I, while Figs. 6 and 7 show the 
agreement that was obtained with the experiment. A 
number of other comparisons have been made, including 
that for a complete cycle of oscillation, and all of them 
show excellent agreement. 

While a check had been established between the measured 
response and conventional theory, a far more significant 
posstbility—that of obtaining the derivatives directly from 
the flight tests—had been opened. Thus, a method for 
determining these aerodynamic constants at high Mach 
Numbers, or under other conditions beyond the scope 
of wind-tunnel data and theory, would be extremely 
worth while. If the equation previously given for la| 
/é is divided by that for |zo| /6, a simple relationship re- 
sults from which Z,, or the slope of the airplane lift 
curve dC,/da, may be obtained directly. The other 
derivatives, including the controlling moment, M;, 
could be extracted from the amplitude and phase rela- 
tionships, provided the test data were of a sufficient 
accuracy. Such a process had been contemplated for 
the present tests, but the instrumentation inaccuracies 
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Fic. 6. Comparison of theory with experimental variation of 
pitchfangle amplitude with frequency of elevator oscillation for 
B-25Jairplane at 175 m.p.h. equivalent air speed. 
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Fic. 7. Comparison of theory with experimental variation of 
normal translational acceleration of the center of gravity with 
frequency of elevator oscillation for the B-25J airplane at 175 
m.p.h, equivalent air speed. 
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Fic. 8. Effect of Z,, My, and M, on theoretical variation of 
normal translational acceleration of the center of gravity with 
frequency of elevator oscillation for the B-25J airplane at 175 
m.p.h. equivalent air speed. 
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Fic. 9. Basic response data with various automatic pilots, B-25J airplane. 


(referred to previously) and resultant lack of phase data 
made this prohibitive. 

In lieu of the above possibility and to establish ac- 
curacy limits for future work, response curves were cal- 
culated in which each derivative was separately set 
equal to zero. The resultant displacement of the re- 
sponse curve gives a practical idea of how closely it 
would have to be established experimentally if the 
derivatives were to be extracted. Fig. 8 shows the 
result of such an investigation on the normal ac- 
celeration response for the derivatives Z,, M,, and 
My. 

The conclusion drawn from studies of this kind was 
that a high but not impractical order of accuracy would 
be required in experimental work directed at this ob- 
jective. 


Concluding Remarks 


While it is desirable to reduce the response measure- 
ments to basic aerodynamic stability and control param- 
eters, direct uses suggest themselves, such as the compari- 
son of aircraft through nondimensional response param- 
eters. Many flying qualities, now measured and speci- 


fied only for quasi-static conditions, stand to be general- 
ized by determination over the complete range of con- 
trol operating frequencies. In fact, sinusoidal testing 
may be a logical way of expressing those items of be- 
havior which vary with pilot rate of control application. 
(From the present tests stick force/g may be plotted 
vs. frequency and extrapolated to zero frequency to 
check the stick force/g as normally measured in wind-up 
turns and steady pull-ups.) It is possible that the con- 
cept of frequency may lead to a more adequate defini- 
tion and understanding of maneuverability. For auto- 
matic control work, simple functions of the measured 
response appear most desirable. 

For conventional aircraft the longitudinal motions of 
primary operational import are those that occur at sensibly 
constant speed. At the lowest frequency of the present 
tests the speed change from trim was less than 2 per 
cent. A fixed control airplane in symmetrical flight 
definitely behaves as a two-degree-of-freedom system, 
whose motions cannot be approximated by the assump- 
tion of one-degree. This was convincingly shown by 


comparing theoretical curves based on both assumptions 
with the experimental data. 


In future tests, the fre- 
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quency range was increased with but slight change in 
speed variations. 

Since the results of these tests are all referred to the 
elevator motion itself, they correspond to conventional 
fixed control stability. Although measurements were 
later made of total normal horizontal tail load (by a 
complex strain gage setup) and elevator hinge moment, 
thus enabling the approximate determination of C,, 
and C;);, the reduction from fixed to free was laborious. 
A more desirable approach to the problem of hinge 
moment response has since been discovered. 


Part III—SuBSEQUENT LONGITUDINAL DEVELOPMENTS 


Introduction 


In the last 1'/. years the longitudinal work outlined 
in the preceding part of this report has been continued. 
In accordance with the original objectives (2 and 3) os- 
cillation tests have been completed with two additional 
automatic pilots—-the Minneapolis-Honeywell C-1A 
and the Sperry A-12—and a considerable fund of longi- 
tudinal response information has been amassed at 
various center of gravity locations, ete. This work has 
all been accomplished on the B-25]J airplane, although 


at the present time a Douglas A-26 airplane is com-_ 


- pletely instrumented and will be used, in part, to obtain 
additional and comparative information. Simultane- 
ously with the testing program, methods have been de- 
veloped for reducing the response data in terms usable 
for automatic control and servomechanism design, 
with continued effort on its utility in the straight aero- 
dynamic design of aircraft (in line with objectives 4, 5, 
6, and 7). 

It will be appreciated that, as an overall endeavor 
of this sort grows in magnitude and as new minds are 
brought to bear on the problems, the number of ideas 
and possible directions for investigation fan out rapidly. 
In the present part of the report an attempt will be 
made to organize and discuss only those ramifications 
that currently promise to fit into the permanent pic- 
ture. These developments were participated in by a 
number of people, and, where practical, direct acknowl- 
edgment is given. 

No attempt will be made to evaluate the relative 
merits, in the present application, of the various auto- 
matic pilots. Such comparative curves as have been 
chosen for this report will not enable any conclusions, in 
this respect, to be drawn. 


Experimental Oscillation Data 


A set of basic longitudinal response curves, showing 
test points, is presented in Figs. 9a—9f. Insofar as the 
data were available, points are shown for the different 
automatic pilots. While in the earlier tests the response 
was measured only in terms of the parameters la! /6 and 
| /é, improvements and additions in instrumentation 
now enable the inclusion of the elevator hinge moment 
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response C,,, and the phase angles between 6 and 6 
(i.e., dg) and elevator hinge moment and 6 (i.e., dyarz). 
It will also be noted that the frequency range is extended 
at both ends from that used in the original tests; other- 
wise the instrumentation setup is similar, though im- 
proved in detail (see Figs. 2), to that used previously, 
The information for the six response parameters is, of 
course, obtained simultaneously for a particular flight 
condition, in this case level flight at 175 m.p.h. equiva- 
lent air speed, 10,000 ft. altitude, with a c.g. location of 
26 per cent M.A.C.—the ‘‘standard”’ condition for the 
B-25 program. 

In the case of the phase angles, ¢» and ¢,, good agree- 
ment was obtained when two totally different types of 
automatic pilot were used, and a number of months had 
elapsed between the tests. This agreement is perhaps 
even more remarkable when consideration is given to 
the size of the filter correction which it was necessary to 
apply to the phase angles. In the case of la] /6 and \n|/6 
the agreement is not so close as could be desired, par- 
ticularly in the low-frequency range. In Fig. 9, curve 
A has been given to illustrate the care that must be 
used in research of this type. From an examination of 
the regularity of the points, one might be led to believe 
in the validity of the data. It was suspected of error, 
however, before the present comparison was made, and 
a Section of the frequency range was repeated with re- 
sults shown at B. The peculiar hump in the C-1A, 
la!/6 response, which initially gave rise to many con- 
jectures, was traced to an oscillator detail. The hinge 
moment phase angle, ¢yzarg, Shows a small “‘lead,’’ which 
is the result of the C,,, component. 

A study of the experimental data in the light of theory 
has been of the greatest assistance in establishing its 
validity and, hence, in assuring a satisfactory test 
technique. It is, of course, generally possible to estab- 
lish an approximation to the parameters for static (o = 
0) and for infinite frequency conditions. 

A second set of longitudinal response information is 
presented in Figs. 10a-10f, where data are shown for 
several air speeds at constant power. The data for 135 
m.p.h. equivalent air speed were therefore obtained in 
climb, while the data for 200 m.p.h. equivalent air 
speed were obtained in dive. These curves were well es- 
tablished by numerous points as for the preceding set. 

The two sets of data have been largely shown to illus- 
trate the amount of information oscillation runs yield 
and the manner in which the basic data are conven- 
iently plotted and faired. Plotting the parameters di- 
rectly against frequency is not, however, always useful 
in judging the effect of an independent variable. In 
Fig. 10b the effect of speed on |n|/5 is indicated satis- 
factorily, but in Fig. 10e the curves are not so simply 
displaced. The use of a frequency ratio as the abscissa 
may be desirable for some cases, such as for data at 
various center of gravity locations. 

To illustrate further the wide range of frequencies 
over which the present tests have been run, the fre- 
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Fic. 10. Basic response data showing effect of varying speed, B-25J airplane. 


quencies for the fixed control short and long (phugoid) 
natural oscillations of the airplane have been indicated 
on Figs. 9 and 10. 

An item of practical importance in obtaining data of 
this type is the reduction of hash in the oscillograph 
galvanometer circuits. Electrical filters are used to re- 
duce this hash and have response characteristics of their 
own which must be known to apply properly corrections 
to the data. This ts illustrative of the necessity of under- 
standing the dynamics of all the systems involved in re- 
search of this type. It appears to be possible to design 
filters whose characteristics in amplitude and phase, in 
the test range of frequencies, are such that corrections 
to the data may be omitted. However, if this is not 
possible, then calibrations and corrections must be made. 
In the case of transient responses, the correct filter de- 
sign is even more paramount, since the corrections are 
difficult to apply if the filter responses are large and 
nonlinear. 


Theoretical Developments 


Once the experimental longitudinal response data have 
been adequately established, the matter of its interpre- 


tation and use arises. Faced with this problem in the 
case of a specific application, a first reaction was to re- 
duce the data, if possible, to the derivatives, whose 
values could then be inserted into the classical equa- 
tions. In the application concerned—a feedback sys- 
tem—it was, however, obviously desirable to express the 
response as simply as possible, and, physically, it was 
felt that for the constant speed motion a simpler means 
must be available. 

E. V. Laitone® has suggested an approach that has 
proved extremely useful in this connection. The solu- 
tions for 6 and n (as given under Analysis and Discus- 
sion, Part II) for motion at constant speed indicated 
that each varies directly with the applied elevator 
moment 61/;. This being the case, it is possible to 
combine the two equations of motion, rather than resort 
to the usual simultaneous solution. The final relation- 
ships for the pitching angular velocity, g, and the result- 
ant normal acceleration, 1, are 


+ bg + kq = — 
ni+bn+kn = UM;Z,6 


where 


b = —(Z, + M, + UM:) 
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k = Z,M, — UM, 
n = ZW 


The quantities in these equations are as previously 
defined. It will be noted that the d and & of the present 
equations correspond to the constants —C, and C; and 
have been resymboled because of their obvious analogy 
to the damping and spring constants of a simple spring- 
mass-resistance system. Since the aerodynamic forces 
and moments are basically assumed proportional to 
velocities, the equations are of one higher order than for 
the corresponding simple mechanical displacement 
system. (Comparable equations have also been de- 
rived for the lateral motion constrained to two degrees 
of freedom—as, for example, by sensing devices and 
automaticcontrol to maintain sideslip, rolling, or yawing 
velocity zero.) 

Once the b and k—which aerodynamically define the 
motion—have been obtained from experimental data, it 1s 
unnecessary for many purposes to accomplish further re- 
duction. They may, for example, be used directly in the 
response equations or for computing the step function re- 
sponse, which in turn enables the evaluation of response 
to arbitrary control motions through Duhamel’s inte- 
gral. The parameters } and k should largely satisfy the 
wants of the automatic control designer for the fixed 
control aspects of his problem. 

The effective damping and spring constants, b and k, 
define the longitudinal fixed control short period of the air- 
plane; thus, 


Prec. = 2n/Vk — (0/2)? 


For the B-25J, center of gravity at 26 per cent, the 
value of b is 3.7 and k = 4.3. Calculations made using 
the B-25J as a standard, and varying design parameters 
such as tail length, horizontal tail area, etc., indicate 
that the dimensional b and k may vary between 2 and 7 
for this class airplane. A 20 per cent increase in tail 
length for a given “‘tail volume”’ would critically damp 
the motion that, for most satisfactory airplanes, is prob- 
ably never far from being critically damped even for el- 
evatorfree. Thespring constant k is seen to be not only 
a function of the static pitching moment slope but also 
of the damping in pitch and the slope of the lift curve. 
If the center of gravity is moved progressively aft, the 
value of k will become zero at approximately 39 per 
cent for the B-25J, which is aft of the fixed control 
neutral point and in the vicinity of the sea level ‘‘ma- 
neuver”’ point. 

For comparison purposes the nondimensional form 
of the equations may be more useful. These are given 
below: 


= bm;(d — 2y) 
= OM 
= + m, + 


(d? + bd +h 
+ bd + 
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where 
b = br = nondimensional, 
k = kr? = nondimensional, k 
n = n(r/U) = nondimensional, 
d = Dr = (d/dt)r 
aw = m/pSc = 1tU/c = relative density parameter 
t = m/pSU = time constant 
ms = 
= 
= TM, 
=\ 
m; = cM, 
¢ = mean aerodynamic chord 
m = mass of airplane 


Other quantities as previously defined. 


Primary Analysis of Experimental Data 


The best means for fairing the experimental response 
data and for determining the values of } and & is the 
“Circle Diagram.’’ This is a graphical scheme origi- 
nally developed in the electrical engineering field and dis- 
cussed by Kennelley.® It was suggested for the present 
use by I. G. Ross. Only the simpler form of this dia- 
gram is discussed here—in relation to the normal ac- 
celeration response, m. A more complex form has been 
used in the analysis of the pitching velocity response, q. 

The variation of normal acceleration response to a 
sinusoidal elevator forcing function is given by 


i + ba + kn = sin ot 


The solution for m, obtained from this equation, may - 


be differentiated and simplified to give: 

UMs 

WV [(k/w) — wf? + 


sin E + tan-! 


b 


where tan-'{ [(k/w) — w]/b} = o; = phase angle be- 
tween 7 and 6. ; 

The maximum amplitude of 7/ \6| for a given fre- 
quency may then be shown to be 


This is the polar equation of a circle passing through 
the origin and having its center on the polar axis, 
where |7|/|6| is the length of the polar vector, ¢, is the 
polar angle, and U\/;Z,,/b is the diameter. It can be 
shown that for ¢; = O (when the vector radius is the 
diameter of the circle), the frequency, w, is equal to the 
undamped natural frequency of the system, wo and k = 
wo”. 

The value of 5 can be found from 

b = (k/w’) — w’ 


where w’ is the circular frequency at which ¢, = +7/4. 
These relationships are illustrated in Fig. 11. 
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Fic. 11. General properties of the circle diagram. 
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Fic. 12a. Circle diagram presentation of experimental data: 


|n|/|6| and gns as a function of w. 
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Circle diagram presentation of experimental data: 
|n|/|Cn| and gnc; as a function of w. 
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Fic. 13. siteseinua variation of stability parameters with 
center of gravity position, B-25J airplane, 175 m.p.h. er 
air speed, 10,000-ft. pressure altitude. (a) Rand b. (b) k’ andjb’. 

In actual use the values of |n|/|5| are obtained di- 
rectly from the experimental data. The theory shows 
that these can be converted to values of |7| /|6| by multi- 
plying by the frequency, w. The phase angles ¢, are 
determined from the experimental values of ¢, by 
adding 90°. Pairs of values of |7|/|é| and ¢,; are then 
plotted, and a circle is faired through the points. Once 
the circle has been established (and it is, of course, of 
great assistance to know that the experimental data 
should fair in this manner), the values of }b and k are 
determined from the simple relationships given above. 

The experimental data when treated in this manner 
show two interesting features. First, the diameter of 
the circle that best fits the data does not fall on the polar 
axis, and the angular difference is in the nature of an 
aerodynamic lead or lag between the elevator angle and 
forcing function moment. The assumptions involved 
in plotting and reducing the data about this displaced 
axis are being further examined. Second, for very low 
and very high values of the frequency, the data departs 
from the circle. In the low-frequency range the de- 
parture may be partly associated with the speed change 
noted in flight; at the high frequencies, where the speed 
is constant, it is probably indicative of the fact that 
the simple velocity derivatives are no longer sufficient 
to describe the motion. Over the large portion of the 
frequency range, wherein the data fair adequately, as- 
surance can be had that: (1) band k exist, (2) that the 
two-degree-of-freedom equations are in line with the 
physical facts, and (3) that the classic velocity deriva- 
tives are probably also existent. 

Fig. 12a is an actual reduction of experimental data 
by the circle diagram method. (Fig. 12a, which is for 
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“fixed’’ control response, should not. be confused 
with Fig. 12b, the hinge moment response, to be 
discussed later.) The angle A is the lag mentioned pre- 
viously. 

Fig. 13a is a plot of experimentally determined ) and 
k values for a wide range of center of gravity locations. 
Each point on this graph represents a reduction by the 
circle diagram; the graph as a whole therefore repre- 
sents a substantial amount of flight testing. It will be 
noted that b is substantially independent of center of 
gravity location, as would be expected. 


Determination of the Aerodynamic Derivatives 


The results obtained from the use of the circle dia- 
gram for 6 and k are indicative of the plausibility of re- 
ducing the data to the familiar aerodynamic derivatives. 
In the frequency range, wherein the circle diagram fits 
the experimental data, such a reduction would logically 
utilize faired values of b and k (vs. center of gravity 
location), which could then be thought of as intermedi- 
ate parameters whose initial fairing increases the accu- 
racy of further reduction. 

Subsequent to the establishment of the parameters 
b and k, math_inaucical analysis yielded similar con- 
stants *’ und k’, which defined the response of the air- 
>" ue to applied elevator hinge moment. These param- 
eters are discussed further under ‘‘Miscellaneous De- 
velopments” but are mentioned here because they are 
not only useful in isolating the aerodynamic derivatives 
but in enabling the determination ‘of other useful aero- 
dynamic parameters. 

The procedure used in the reduction will now be de- 
scribed. (For convenience in handling the data, all re- 
sults are in terms of the dimensional quantities.) 

It has been previously noted (Part II) that a simple 
relationship for Z, may be found by dividing the @ and 
w responses. In terms of the more convenient param- 
eters n and q, this relationship is: 


= U/V1 + (@/Zy)? 


Fig. 14 shows Z,, as a function of frequency, from ex- 
perimental data at a number of center of gravity loca- 
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Fic, 14. Variation of the derivative Z, with frequency of 
oscillation, B-25] airplane, 175 m.p.h. equivalent air speed, 
10,000-ft. pressure altitude. 


TABLE TI 


COMPARISON OF ESTIMATED VALUES OF STA- 
BILITY DERIVATIVES AND PARAMETERS WITH 
VALUES CALCULATED FROM FLIGHT TEST DATA 
6-25J AIRPLANE , |75MPH.EAS. , 10,000 FOOT ALTITUDE 


C.G. POSITION 6.6. POSITION 
26% MAC | 28% MAC | 26% MAC | 28% MAC 
4Zu | - 1.00 * | -1.00 -1.25 -1.08 
My ‘Sec. |-0.0156 * | -0.0125 
Ms Yer -0.0023 | —0.0023 | —0.0015 —0.0007 
‘ec -1.55 -1.55 —1.60 -2.20 
M; Mec? | -8.64 | —8.64,%, | —8.00 
F- — 0.45 0.45 0.75 0.70 
Cys Yoec. |-0.0039 |-0.0039 » | -0.0054 | —0.0054 
Choc [-0.0035 |-00035 * | -0.0046 | ~0.0041 
k 6.21 5.28 6.90 6.10 
b ‘Kec 3.24 3.24 3.50 2.50 
ki ‘Lec? 1.27 0.34 4.80 3.70 
b' Kec. 2.27 2.27 260 | 2.60 


* OBTAINED FROM WIND TUNNEL TEST DATA. 


#®*s OBTAINED FROM STATIC STABILITY FLIGHT TEST 
MEASURE MENTS. 


NOTE: SOME VALUES GIVEN IN THIS TABLE WILL NOT COMPARE 
EXACTLY WITH THE CORRESPONDING VALUES OF TABLE TI. 
IMPROVED DATA WERE AVAILABLE FOR THE CALCULATIONS 
SUMMARIZED IN THIS TABLE 


tions. These data are again for the B-25J airplane. 
From an examination of this plot, —1.25 was taken 
as the most probable value for 26 per cent c.g. in the 
low to moderate frequency range. 

The value of 1/; may now be determined from the 
diameter of the fixed control circle diagram which 
equals UM;Z,/b. The corresponding diameter of 
the hinge moment circle diagram (the |7!/C, circle) 
is equal to UM;Z,,/C,,’ and permits of the isolation of 
Crs- 

Subtraction of the mathematical relations for 6’ and 
band k’ and k yields two equations that may be solved 


de 
simultaneously to give C,, and da’ thus, 
a 
(b’ — b) + (l/U)(R’ — k) 
Chg U U 


da 


where / = tail length (center of gravity to aerodynamic 
center of horizontal tail) and M; is per degree. 

The damping in pitch is found from the definition of 
k if M, is known from wind-tunnel or static flight data, 
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while 7; results from the definition of 0; thus, 
ky?/\OC, 
My, = —(6+ + M,)/U 
where k, = radius of gyration about y axis. 

It will be noted that a value of M, is required from 
other sources for the complete reduction of the oscilla- 
tion data. Further developments may make even this 
unnecessary. 

The results of an analysis (for 175 m.p.h. equivalent 
air speed and a center of gravity location of 26 per 
cent) are shown in Table II, where comparison is made 
with values estimated from theory and powered model 
tests. A reduction has also been made from data at 
28 per cent and is included. An examination of this 
table shows that the test results are, in general, reason- 
able. 

The significance of the above results lies in the demon- 
stration of the feasibility of the new method in isolating 
aerodynamic derivatives, and other basic aerodynamic 
data. Some of these data are unobtainable from flight 
by other means, and the overall method and procedure 
for obtaining the derivatives is as simple as any that 
has yet been proposed. In fact; once 4, k, etc., are 
known, the reduction to the aerodynamic parameters 
may be completed in less than an hour. 

Obviously, much remains to be done before the 
method can be considered for routine engineering work. 
In the first place, a thorough examination of the reduc- 
tion accuracies is being made to ensure that the present 
method is optimum. In the second place, it is proposed 
to reduce a large amount of flight-test data by the 
method, thus enabling checks and accumulation of ex- 
perience. 


The Step Function 


In addition to the sinusoidal, there are many other 
forcing functions that may be of practical importance in 
stability and control work. Thus, for many years it has 
been customary to examine, for certain purposes, the 
result of abrupt manual control deflection. A further 


study of this technique, which will be referred to as the 
“step function,’ has been made and will now be dis- 
cussed. 

First, it should be emphasized that the sinusoidal 
forcing function has the great advantage of yielding a 
large volume of steady-state data, and it is likely that 
it will always excel other types of forcing functions in 
this respect. On the other hand, it has real disad- 
vantages in the flight time required and the problems of 
data reduction— both of which may prove excessive for 
routine use. Also, there are many rapidly changing 
flight conditions (such as dive, or a maneuver at high 
acceleration) wherein time is not available for a steady- 
state examination, even if the human pilot would re- 
linquish the flight controls completely and allow them 
to be operated automatically. (It is true, of course, 
that automatic control will be used more and more in 
these critical conditions of flight, but for the present 
they are approached with manually operated controls.) 
The forcing function of short duration, with reduction of 
the necessary data from the transient response, therefore 
looms as an attractive possibility. 

A proposed technique will be discussed briefly. Thus, 
suppose that a pilot could, by pushing a button (lo- 
cated on the control stick), momentarily lock the cock- 
pit control, while at the same time the control surface 
was displaced through a small angle (0.25° to 0.5°) rela- 
tive to the cockpit control. Records of the airplane 
response in this brief period could be simultaneously ob- 
tained, and, if the information were sufficient and satis- 
factory, the imagined technique would be exceedingly 
useful in rapidly obtaining data in all sorts of maneuver- 
ing (and steady) flight conditions. Hydraulic and pneu- 
matic mechanisms can easily be visualized which would 
make such a technique workable. In comparison to 
the sinusoidal method, it would be possible to simplify 
the instrumentation and eliminate the automatic pilot, 
which would be desirable for fighters and other 
“cramped quarters”’ airplanes. The overall time for the 
step application could be under the control of the pilot 
or adjusted to some predetermined amount (on the 
order of 0.5 to 0.8 sec.). 


The establishment of the technique described above depends primarily upon its relationship to the mathemat- 
ical step function, which assumes that the control is displaced instantaneously and that the corresponding ap- 
plied aerodynamic moment acts in the same manner. By straightforward analysis the response to a unit elevator 
step (6) 1(t), in terms of the dimensional 6 and k, may be expressed as 


n(t) _ 
g(t) | + (2k/Z..)}? V4k — 


where, b? < 4k (oscillatory) 
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Fic. 15. Comparison of measured and calculated response to 
a step deflection of the elevator, B-25J airplane, 175 m.p.h. 
equivalent air speed, 10,000-ft. pressure altitude, center of 
gravity = 26 per cent M.A.C. 


There was no assurance whatsoever that an equiva- 
lent of this mathematical step could be obtained in 
flight or even approached. Nevertheless, from a prac- 
tical standpoint the possibility that an experimental 
step function could be used directly with the Duhamel 
integral to compute responses to arbitrary control mo- 
tions (and that for certain cases good, if not better, 
agreement might be obtained than from a calculated 
step, based on experimental b and k) could not be over- 
looked. Obviously, a number of comparative checks 
would have to be made; thus: 

(1) A comparison of a step function calculated from 
sinusoidal b and k data, with an experimental step. 

(2) Acalculation of b and k values directly from an 
experimental step, and a comparison with b and k 
from experimental oscillation data. 

(3) A comparison with the results of an experimen- 
tal arbitrary maneuver by (a) the use of b and k from 
sinusoidal tests and the Duhamel integral, and (b) the 
use of an experimental step and the Duhamel integral. 
If the experimental arbitrary maneuver results could 
be reproduced by either (a) or (b), then the validity of 
the particular technique has been established. 

Experimental steps have been performed using both 
manual and automatic control. In most cases the 
elevator was snapped up about 1° to 1!/2° from trim, as 


' fast as possible manually, or at the maximum servo- 


output rate. In each case the control was displaced be- 
fore any change in attitude of the airplane was recorded, 
although an extremely small momentary downward 
acceleration from 6Z; was noted. From such tests as 
were made with other control surface displacements, it 
appears that for small angles the response is linear with 
displacement. From practical flight aspects the tech- 
nique was entirely satisfactory. 

A comparison as suggested under (1) is made in Fig. 


‘15. The caiculated curve was obtained by inserting 


measured 6 and k values in the step function response 
equation, given above. The experimental curve has 
been modified only by the application of electrical filter 
corrections, which in this case required a Fourier analy- 


‘ 


sis. At the end of 0.8 sec., the airplane attitude had 
changed appreciably, and, since the speed had started 
to drop, the comparison is no longer entirely valid above 
this time. (With a “push down’ it appears possible 
to reach an asymptotic value before appreciable speed 
change.) If M;, is not exactly in phase with 4, as is as- 
sumed in the equations, and if, in addition, the b and 
k values from the sinusoidal reduction do not ade- 
quately take into account the effect of 6Z;, a shift of the 
time scale of the measured curve may be required be- 
fore a comparison is made. If the measured curve is 
arbitrarily shifted to the left by, in this case, approxi- 
mately 0.1 sec., the agreement up to 0.8 sec. is then re- 
markably good, the general shapes comparing well and 
the slopes being identical in the range 0.4 to 0.5 sec. 

A comparison of the time and effort required to ob- 
tain the two curves is interesting. The sinusoidal data 
for b and k took about 3 hours of flight, while securing 
the step data was a matter of less than5min. The time 
for record processing and data reduction is, of course, 
much greater for the sinusoidal case and is of the order 
of ten to one. 

The step function results to date are inconclusive. 
Further tests are contemplated with a series of dif- 
ferent step magnitudes (using improved filters) to 
establish empirically the initial shape of the curve. 
Meanwhile, methods of obtaining 6 and k from a step 
function response are being examined. The asymptotic 
value determines kif W; and Z, are known; the extrac- 
tion of b is more complicated. The practical aspects of 
the entire matter rest upon the arbitrary maneuver, 
comparison (3), which has recently been performed but’ 
as yet not analyzed. 

In general, control characteristics (particularly hinge 
moment) are a function of the rate of application, and it 
is now well known that adverse effects may be obtained 
during rapid maneuvers. The automatic “‘step func- 
tion equipment” enables a direct study of these effects 
by obtaining data at a number of predetermined con- 
trol rates—i.e., this equipment may be used to apply 
the control at variable rates if desired. 


Miscellaneous Developments 


Speed Changes 


Some further discussion of the constancy of forward 
speed during the oscillation tests is required. In Figs. 
16a and 16b are shown measured longitudinal accelera- 
tions (properly reduced by removing the pitch angle 
component of gravity) and the indicated variations in 
speed from the air-speed indicator. (The latter are 
merely indicative of what the flight crew could ob- 
serve directly.) Although the acceleration data are con- 
clusive, extensive calculations have been made at both 
Massachusetts Institute of Technology’ and Cornell to 
examine theoretically the probable speed variation over 
the frequency range. If values of b and k are computed 
from expressions that include all the U velocity terms 
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and are compared to the computed two-degree-of-free- 
dom case, the following results are obtained (B-25]J air- 
plane, 175 m.p.h. equivalent air speed) : 


b k 
Two degrees of freedom 3.78 4.284 
Three degrees of freedom 3.77 4.252 


Speed change effects may be taken as less than 2 per 
cent at frequencies as low as 0.5. 

In connection with the effect of speed changes, etc., 
it should be noted that, for the cases where the experi- 
mental points fall off the circle diagram (at high or low 
frequencies), an empirical method has been developed 
for obtaining the effective b and k values, which are, of 
course, functions of frequency. 


Hinge Moment Response 


The hinge moment response may be expressed in rela- 
tionships comparable to those for the control angle 
response, and effective damping and spring constants b’ 
and k’ may befound. A circle diagram of this reduction 
for a particular case is given in Fig. 12b, and the experi- 
mental variation of these parameters with center of 
gravity location is shown in Fig. 13b. It should be 
noted that the hinge moment data are obtained from 
strain gages at the control horns. Since servomotors in 
automatic pilots have limited torque outputs, the auto- 
matic pilot manufacturer is interested not only in the 
response of the airplane to control displacement but 
also to applied hinge moment; 0b’ and k’ fill this latter 
need. The automatic pilot may be thought of as a de- 
vice that is neither ‘‘fixed’’ control nor “‘free’’ control, 
and the aerodynamic specifications for the automatic 
pilot must then include’, k, b’, and k’. 

From the aerodynamicist’s standpoint, the hinge 
moment parameters offer additional relationships that 
can be used, as in the case of the derivative analysis, to 
isolate design factors of interest. 


Special Sensing Equipment 


In the course of the program a number of special 
sensing units have been developed or used. Among 
these are angle of attack, angle of yaw, indicated air 
speed, pressure altitude, and gsensing. The output from 
these sensing units is fed into the automatic pilot to 
maintain a given value of the variable sensed or to 
change (or rate of change) it at a predetermined rate. 
A combination of the g and speed sensing enables, in a 
single automatic maneuver, the procurement of all of 
the information desired for SF/g flying qualities tests, 
which has formerly been obtained (and only with consid- 
erable difficulty) in wind-up turns. A slow rate of 


change of control position device (for moving controls 
continuously and as slowly as desired) is of consider- 
able use in obtaining static data from flight tests. 
Curves of static elevator angle, and stick force vs. speed, 
were obtained when this device was set to operate at the 
rate of 1° of elevator per 5 min. Further tests will indi- 
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Fic. 16. Variation of speed during sinusoidal oscillations at 


several frequencies, B-25] airplane, 175 m.p.h. equivalent air 
speed, 10,000-ft. pressure altitude. (a) Calculated from longi- 
tudinal acceleration records measured in flight. (b) Observed 
from photo-observer records of air-speed indicator taken in 
flight. 


cate how fast a control may be moved before the as- 
sumption of static conditions is violated. 


Effect of Magnitude of Control Motion During 


Oscillation Tests 


It was initially assumed in oscillation testing that the 
response was linear with control deflection, provided 
the deflections were small. Since it is desirable in some 
respects to use as large deflections as possible, experi- 
mental-measurements have been made to ascertain the 
effect of changing this parameter. No effect on the 
response per unit of control deflection (up to 3.0° total 
deflection, the highest tried) was found. 


Static Flight Checks 


A number of static parameters (Cz,, Cng, etc.) have 
been measured by well-known methods to check the in- 
formation from oscillation tests. This type of checking 
has been a continuous part of the program and has been 
in part referred to in the derivative analysis. 


Notes on the Instrumentation 


Unlike many types of flight testing, a large part of 
the instrumentation for the present program has been 
especially developed. This instrumentation is largely 
electronic, and, while it is not the purpose of this re- 
port to go into detail, a brief discussion of some of the 
problems is given below. 

Perhaps the major instrumentation problem to date’ 
has been that of filters. The necessity for adequate fil- 
ters has been shown by a sample of unfiltered oscillo- 
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Slow rate of change of control position mechanism 
used with A-12 autopilot, B-25J airplane, 


Fic. 17a. 
graph record in Fig. 2, and the necessity for an under- 
standing or elimination of the filter responses has been 
discussed. 

During the course of the tests, a phenomenon ap- 
peared in the records which has been referred to as the 
““‘beat.’’ The oscillograph records showed, and it was 
visually confirmed in flight, that the airplane did not 
oscillate at regular amplitude but at even intervals 
made a cycle of higher amplitude. Although every ef- 
fort was made to isolate the cause of this phenomenon, 
it suddenly disappeared and has never re-occurred. 
During its life, however, special instrumentation and 
modification of the oscillator and automatic pilot were 
required for its study, which posed original problems for 
the Instrument Section. 

Fig. 17a shows the device developed for moving the 
elevator at extremely slow predetermined rates. This 
is an ingenious adaptation of an ordinary alarm clock. 
Fig. 17b illustrates the equipment for acceleration 
sensing and for applying step functions. The: design 


normal acceleration sensing equipment. 


Fic. 17b. Top: 
Both used with A-12 auto- 


Bottom: step function equipment. 
pilot, B-25J airplane. 
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requirements of the former were naturally stringent to 
avoid any possibility of exceeding allowable accelera- 
tion limits and, hence, have structural failure of the air- 
plane. The device incorporates a unique g limiting 
device that can be present to a safe allowable value. 
With appropriate servos, the step function equipment 
can be used to obtain step functions of any magnitude. 
In operation, the control button is pushed when the 
necessary flight conditions have been established. The 
control snaps over and remains deflected until the but- 
ton is released. 

Instrumentation development, calibration, and main- 
tenance for a research of this type must be looked upon as 
a major part of the operation and the program. 


Weight, Center of Gravity, and Inertia Control 


An important practical consideration in this type of 
flight testing is that of the control of weight and its dis- 
tribution. Early in the program, water transfer pro- 
visions were made so that the c.g. could be shifted in 
flight. A tail water tank (400-lb. capacity) allows for a 
shift of 5 per cent M.A.C. at normal gross weight. The 
water in this tank can be controlled to +10 lbs., and 
the c.g. may be placed in flight to a precision of 0.2 per 
cent M.A.C. The weight at a particular time may be 
ascertained to about 1 per cent. 

Basic moment of inertia data as supplied by the 
manufacturer had been used with appropriate correc- 
tions for configuration changes. Calculations revealed 
that large moment of inertia changes (on the order of 
30 per cent for the pitch axis) were encountered on 
shifting water ballast for c.g. control. Under these 
circumstances, in consideration of the dependence of 
the response measurements on’ moment of inertia, a 
study was made of methods of experimentally deter- 
mining the moments of inertia. Jt was concluded that the 
method of swinging, which has occasionally been used, 
suffers inaccuracies unless care 1s taken to use the optimum 
pendulum length for each case. The practical problems 
connected with locating a supporting structure of the 
correct height and swinging an airplane of the size of 
the B-25 led to the development of a ground oscillation 
technique. The airplane was pivoted at the wing-jack 
points, while supported at the nose gear location on a set 
of coiled springs of known characteristics. Records 
were taken on the oscillograph of the time history of 
the pitch angle during a free oscillation. For practical 
purposes it was found possible to neglect the damping 
in the system, since several cycles of oscillation would 
be obtained with little change in amplitude. This 
facilitated the determination of the longitudinal mo- 
ment of inertia from the data. The complete results of 
these tests were not available in time for inclusion in the 
present report. 

Consideration is also being given to the development 
of techniques for the experimental determination of the 
moments of inertia about the X and Z axes of the air- 
plane for the lateral program. 
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Part IV-—FuTuRE DEVELOPMENTS 


Dynamic Lateral Stability and Control 


A program of dynamic lateral exploratory tests has 
been completed on the B-25 airplane. These tests, 
which utilized automatic control, involved about 20 
hours of flying. In addition, extensive plans, theoretical 
developments, and the preparation of special instrumen- 
tation have been under way for the complete lateral 
program. 

In general, the philosophy of the lateral program has 
been that, once the longitudinal motion was thoroughly 
mastered, the technique and experience could be trans- 
ferred to the lateral case. In this case, involving as it 
does the unsymmetrical motion of the airplane, auto- 
matic control offers possibilities in constraining a par- 
ticular velocity component, resulting in an effectively 
simpler motion for study and analysis. Thus it is ap- 
parent that the sideslip velocity might be maintained at 
zero through the use of a sensing device (controlling 
the rudder through the automatic-control) based on yaw 
angle and/or lateral acceleration and, if necessary, their 
derivatives. Such constraint would reduce the motion 
to that of a ‘“‘two-controlled’”’ type (such as would be 
experienced with an airplane of high static directional 
stability and damping), since the sideslip velocity 
would be maintained approximately zero at all times. 
Many possible types of constraint, of which that de- 
scribed is but one, may be visualized for the lateral mo- 
tion by the appropriate adaptation of automatic con- 
trol. 

The exploratory tests have also included examination 
of a wide variety of suggested oscillation techniques, 
with, in each case, the primary objective of establishing 
their feasibility from the practical flight operating 
standpoint. While space does not permit a complete 
review of these many tests, it can be mentioned that 
provision was made for oscillating the rudder and ai- 
leron separately or in combination. These control 
oscillations could be performed with various functions 
of the automatic pilot in effect. Thus, the control 
oscillations could be superimposed upon the normal 
automatie pilot function (which is that of maintaining 
an overall trim attitude), as in the longitudinal case, or, 
at the other extreme, the servomotors could be made to 
oscillate the controls with no normal automatic pilot 
function present. Further provisions were made for 
utilizing various amounts of gyro signal (adjustable 
in flight) in the rudder and/or aileron circuits. The re- 
cording instrumentation was maintained at a minimum 
during the exploratory tests, although a small photo- 
observer was utilized for the specific purpose of es- 
tablishing that the aileron and rudder motions were 
satisfactorily sinusoidal during the oscillations. The 
various test conditions were set up largely through vis- 
ual instrumentation of sideslip angle, side acceleration, 


etc. 


The theoretical work performed simultaneously with 
the exploratory flights was directed at an examination 
of the various techniques under test, with an eye to es- 
tablishing the lateral information that each would 
yield, the specific recordings that would be required, 
and the methods of data reduction and analysis. 


As a result of all of this exploratory work, the scheme 
that offers the most promise, both from a flight and 
analysis standpoint, may be described as follows: 


Provision is made to oscillate both tke ailerons and 
the rudder simultaneously and at identical frequencies 
by means of the automatic controlling system. In the 
case of either the rudder or aileron control, the normal 
automatic pilot function (gyro signal) is introduced in 
an amount necessary to maintain average trim attitude 
during the oscillation. The control motions are limited 
to small amplitudes, which allow the ensuing lateral 
motions to be considered as independent of the longi- 
tudinal motions. (The elevator remains on normal 
automatic control throughout.) Provision is made for 
the manual adjustment of the phase between the two 
control motions and their maximum amplitude ratio. 
By an appropriate adjustment of the phase and ampli- 
tude ratio control, it is possible to show tleoretically, 
and has also been proved experimentally, that three 
special lateral maneuvers may be established—each 
maneuver being defined by the elimination of one of 
the lateral velocity variables, that is, sideslip velocity, 
rolling velocity, or yawing velocity. Thus, the first 
maneuver, in which the sideslip velocity v = 0, consists 
of a series of perfectly coordinated turns; the second 
maneuver, in which the rolling velocity = 0, consists 
of a yawing oscillation in which the wings remain level; 
while the third maneuver, in which the yawing velocity 
r = 0, consists of a rolling oscillation at constant head- 
ing. Each of-the maneuvers is repeated, as in the longi- 
tudinal case, for the entire usable range of frequencies 
(0 to 8 rad. per sec.), and an oscillographic recording is 
made of all of the pertinent variables of the motion, 
such as the angular displacements, the angular veloci- 
ties, linear and angular accelerations, etc., as well as 
the phase angles of the ‘‘forcing’’ control motions and 
the amplitude ratios. 


The theoretical analysis indicates that, from the data 
obtained through tests of the above type, it should be pos- 
sible to reduce, in pairs, all of the lateral stability and 
control derivatives—both static and dynamic. It may even 
be possible to obtain all of this information from two of 
the maneuvers, or it may be possible to obtain duplicate 
information on vaiious derivatives if all of the special 
maneuvers are used, thus enabling an extremely valu- 
able experimental check. The reduction of the data 
can proceed upon semigraphical lines, after the manner 
of the longitudinal analysis and, hence, enable the de- 
termination of effective spring and damping constants 
of those modes of the lateral motion which are of prac- 
tical import. . 
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This ‘‘double oscillation’’ lateral forcing function 
technique, which represents such a powerful attack on 
the lateral case, was first suggested by J. B. Rea while 
fulfilling the requirements for his Doctorate at the In- 
strumentation Laboratory, M.I.T. The development 
of the flight technique is the work of the Flight Re- 
search Department at Cornell Aeronautical Laboratory. 
Advances in the method for extracting the individual 
lateral derivatives from the test data are the result of 
the joint efforts of I. G. Ross, W. O. Breuhaus, and D. 
W. Whitcomb, of the Flight Research staff. 

To complete the lateral picture, it is likely that the 
step function technique will be utilized, and, in addi- 
tion, static flight checks will be performed by stand- 
ard ‘‘applied moment”’ methods. 

In judging the significance of a flight method for the ex- 
traction of all of the lateral derivatives, attention should be 
directed to the fact that at the present time there is no com- 
monly used method, either wind tunnel or flight, for ob- 
taining the damping and coupling terms. Furthermore, 
there is no flight method available for obtaining directly 
the magnitudes of the control effectivenesses and the im- 
portant static lateral derivatives, such as the effective di- 
hedral and the directional stability. 


Combined Longitudinal and Lateral Motion 


Flight measurements of the combined dynamic longi- 
tudinal and lateral motion are almost nonexistent, al- 
though a considerable fund of static information on such 
items as the pitching moment due to sideslip is avail- 
able. When automatic control is used for rapid ma- 
neuvering, such as tracking, considerable interest is at- 
tached to the entry and recovery from turns and the 
characteristics in steady turns with and without side- 
slip. The longitudinal characteristics, thought of 
alone, may be applicably modified by the steady incre- 
ment of normal acceleration, which maintains in steep 
turns or rapid pull-ups. Furthermore, the motions re- 
sulting from rough air are a combination of the longi- 
tudinal and lateral modes, even though the airplane is 
attempting to proceed along a steady rectilinear flight 
path. 

Exploratory tests have been completed which demon- 
strate the practical feasibility of performing elevator 
oscillation tests under accelerated conditions in turns of 
various degrees of bank. Other tests have been run 
which show that similar tests may be made with vari- 
ous degrees of initial yaw. Finally, by a combination 
of the oscillator control and automatic pilot turn con- 
trol, oscillation data are available during turn entry. 
To survey and isolate the critical conditions for com- 
bined longitudinal and lateral motion will be a major 
research effort in itself. 


Long-Range Dynamic Research Program 


A program for the continuance of the research in 
basic dynamic stability and control (which is in ac- 


cordance with the initial objectives of Part I) is mapped 
out in preliminary form below. Such a program would 
represent in the neighborhood of 5 years of research 
and in the author’s opinion is a logical one for furthering 
the cause of dynamic analysis so that its import will be 
effectively felt in several allied fields of design. The 
work outlined in Phase I is now largely planned or in 
work under the sponsorship of the Air Materiel Com- 
mand. In Phase II are collected those items that 
should follow immediately upon the completion of 
Phase I. This phase is largely concerned with the ap- 
plication and “‘pilot’’ trial of methods developed in 
Phase I and is a necessary and important phase in any 
broad research. Phase III is concerned with specific 
and currently foreseeable dynamic problems of suf- 
ficient complexity and importance to warrant particu- 
lar research. It can be anticipated that in this phase 
all the tools of previous research would be effectively 
utilized. 


Phase I 


(1) Perfection of longitudinal methods with particu- 
lar attention to step function. 

(2) Adaptation of methods to lateral and com- 
bined longitudinal and lateral motion. 

(3) Development of sideslip, roll, etc., sensing 
equipment as required by (2) above. 

(4) Study of application of the longitudinal oscilla- 
tion technique to high-speed (subsonic) airplanes and 
to sub-, trans-, and supersonic missiles. 

(5) Study of application of oscillation data to spe- 
cific problems in automatic control and simulator de- 
sign. 

(6) Study of application of advanced electrical engi- 
neering methods to aerodynamic stability and control 
problems. 


Phase II 


(1) Longitudinal oscillation tests on a series of 
widely different types of airplanes for the purpose of 
amassing experimental information over a wide range 
of design variables. 

(2) Simplification of longitudinal oscillation tech- 
nique for routine use and for (1) above. 

(3) Repeat (1), using best lateral technique. 

(4) Study of information from (1) and (3) above 
in reference to (a) the specification of dynamic flying 
qualities characteristics and (b) the specification of com- 
mercial and military airplane behavior on automatic 
control. 

(5) Application of the longitudinal oscillation 
method to a high subsonic speed airplane (P-SO class) 
on the basis of Phase I (4). 

(6) Longitudinal oscillation tests on supersonic 
missiles with telemetered information. 

(7) Extension of longitudinal oscillation technique 
into the unstable range (on a B-25J class airplane). 
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(8) Development of electronic simulator for utiliz- 
ing 6 and k for routine use in the research. 
(9) Development of dynamic design charts for the 
new parameters, k and 3), and in nondimensional form. 
(10) Study of rough-air behavior of airplanes on 
automatic control. 


Phase III 


(1) Research on practical rough-air control equip- 
ment. 

(2) Stall research using automatic control and the 
development of stall control. 

(3) Special missile applications. 

(4) Research in the measurement of the dynamic 
characteristics of helicopters. 


GENERAL CONCLUSIONS 


Substantial progress has been made in the attain- 
ment of the objectives laid down for this research: 

(1) New methods and techniques have been de- 
veloped which have wide application in the measure- 
ment of dynamic stability and control in flight. 

(2) By the use of these methods, a large amount of 
full-scale dynamic information has been obtained on 
service-type aircraft for the longitudinal motion. 

(3) Substantial progress has been made in simpli- 
fying the analytical approach, particularly in reference 
to aerodynamic data for automatic control design and 
in enabling rational reduction of the flight informa- 
tion. 

(4) Some insight has been obtained into the range 
of validity of the classical assumptions, and an ap- 
proximate method for determining the aerodynamic 
derivatives from flight has been demonstrated. 

(5) The use of automatic control for this type of 
work has been thoroughly established, as well as its at- 
tributes in other types of flight testing. 

(6) A tremendous amount of experience in all 
phases of this type of research is now available to be 
directed at the rapid attainment of the remaining ob- 
jectives. 

Note: More detailed discussion of both the theoretical 
and experimental aspects of this research has recently 
been published in a series of Cornell Aeronautical Labora- 
tory technical reports.*—*! For a résumé of lateral re- 
sponse theory and advances in the longitudinal treatment, 
reference is made to these reports. 
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APPENDIX—SYNOPSIS OF PREVIOUS DYNAMIC RESEARCH 


Up to about 3 years ago, at the time of the initiation of the 
present work, progress in dynamic stability and control in the 
English-speaking countries was largely associated with conven- 
tional airplane development. Although in Germany and, more 
recently, in other countries emphasis has been placed on dynam- 
ics in missile and pilotless aircraft application, it is currently dif- 
ficult to evaluate these efforts from the standpoint of fundamental 
progress. A completely digested picture of the German work is 
not yet available, and a general review is still hampered by exist- 
ing classification. A large number of references to basic German 
dynamic investigations, none of them translated and many not 
readily available, have come to the author’s attention. A large 
portion of these relate, interestingly enough, to combined aero- 
dynamic and automatic control problems. The hopelessness of 
evaluating German progress in this field at the present time is ap- 
parent. Some notes on interrogations arranged through Foreign 
Information, T-2, Air Matériel Command, are, however, in- 
cluded at the end of this Appendix. 

A sufficient number of key or typical references are given below 
to enable the reader to pursue a particular phase of the literature. 
No attempt is made at.a complete bibliography. 


Classical Theory and Early Model Experiment 


Bryan’s pioneering work in aircraft dynamics, based on the 
rigid dynamics of Routh, was initiated in the year of the air- 
plane’s birth (1903), and its complete presentation in Stability in 
Aviation’ in 1911 marks the beginning of the so-called. classical 
stability approach. The complete mathematical statement of the 
fixed control motion of an aircraft was presented in six simultane- 
ous linear differential equations with constant coefficients, predi- 
cated on the assumption of small disturbances from steady sym- 
metrical flight and with component air forces as linear functions 
of the velocities (i.e., the air force derivatives). These linearized 
equations are written for axes that are fixed in the airplane. For 
the reason of airplane symmetry and small disturbances, the 
equations normally separate into two groups—the longitudinal 
and the lateral. It is these equations, without fundamental 
change, that have for the last 35 years represented the most 
common approach to the study of airplane stability. Almost the 
sum total of the mathematical and experimental work in aircraft 
dynamics has therefore evolved about equations whose rigorous 
application does not extend beyond small departures from steady 
trimmed symmetrical flight and which are incapable of handling 
problems of finite maneuvering in curved flight—a limitation, 
however, nowhere near so drastic as it might first appear. The 
earliest, and one of the few, exceptions to Bryan’s approach was 
that of Lanchester,? which depends on great simplifying assump- 
tions, largely prohibiting its use in research or design. 

In the years following Bryan’s original presentation, to about 
1920, much effort was expended in further mathematical refine- 
ments and generalization of this theory. In Great Britain, con- 
tributions were made by Bairstow, Jones, Glauert, and Brodet- 
sky,’°- 13 while in this country Wilson,'* in some of the earliest 
N.A.C.A. work, studied the motion of an aircraft encountering 
gusts, pointing out among other things the significance of the 
short period to acceleration response. By 1919 Glauert had pub- 
lished a complete theoretical treatment of the free-control longi- 
tudinal stability and the stability derivatives of the propeller,’® 
while Bateman had introduced a form of the stability diagram. 

In parallel with the theoretical was the comprehensive develop- 
ment of wind-tunnel methods and techniques for the measure- 
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ment of the individual stability derivatives. While this work 
originated in Great Britain, it was introduced into this country 
in about 1915 by Hunsaker, who, with others, ran complete dy- 
namic tests on two models in the old 4-ft. rectangular tunnel at 
the Massachusetts Institute of Technology. These results are 
presented in a Smithsonian Report, Dynamic Stability of Aero- 
planes,” perhaps the best known early presentation of the sub- 
ject in this country and one that can hardly be surpassed asa 
basic reference today. The work was continued by Warner, 
Klemin, and others with increased emphasis on the effect of de- 
sign changes. In Great Britain, the refinement of model methods 
for the determination of the individual derivatives continued, 
though at a decreasing pace, throughout the 1920’s. Glauert, 
Bryant, Irving, Relf, Bradfield, Halliday, and others*~?! ac- 
complished an amazing amount of research with continuous 
rotation balances and simple (one degree of freedom) oscillation 
techniques in the wind tunnel and with the whirling arm. Only a 
small typical part of this work is cited in the references. 

The introduction of a nondimensional form of the stability 
equations by Glauert”* in 1927 was a first step toward their engi- 
neering use, and it facilitated the examination of the effect of pri- 
mary variables upon the motion. Based on these equations, 
Gates?’ surveyed the unstalled longitudinal motion, presenting 
the results graphically in stability diagrams. It is doubtful, 
however, that even this simplified presentation was greatly used in 
design. 

Complete texts of the progress of dynamic stability and con- 
trol were available for the student in 1920 in a paper by Glauert*! 
and in Bairstow,** Applied Aerodynamics, and Wilson, Aero- 
nautics, which were followed by Warner,”” Airplane Design, in 
1927. A review of progress at this stage indicates that, while 
the Bryan classical approach had become well established in re- 
search and model methods were available for experimental de- 
termination of the derivatives, its use in practical design was neg- 
ligible. In fact, the static approach, utilizing only the static 
derivatives, had become the recognized design tool. Some of the 
shortcomings of the dynamic approach, such as its almost com- 
plete confinement to the determination of the period and damp- 
ing of the phugoid and the prohibitive time for a complete wind- 
tunnel test of a particular model, were well recognized. In later 
years, even in research, the specialized balance and whirling arm 
equipment became less and less used, for it was not adaptable to 
the determination of information at higher Reynolds Numbers. 
Other seemingly more profitable avenues of research led to a 
change in emphasis as the picture developed. 


Initial Flight Measurements and Stall Studies 


Almost as soon as wind-tunnel methods were proposed for the 
determination of the derivatives, attention was directed to their 
measurement in flight. Preliminary dynamic flight tests were 
performed by the British prior to 1920,¥% but the outstanding 
work of the early period is that of Norton, Warner, and Allen,” 33 
at the N.A.C.A., which took place between 1919 and 1923. It is 
interesting to note that this work marks the beginning of N.A.- 
C.A. flight-test and instrumentation development and E. T. 
Allen’s debut as an engineering test pilot. In the course of this 
program, a number of the derivatives were determined(l/,, /,, 
Nes Vor My, etc.), and use was made of the ‘‘step function’”’ (sand 
dropped from wing tip) and the ‘known applied moment” 
(parachute at wing tip) techniques. Consideration was given to 
desirable flying qualities in relationship to basic stability and 
control, and comprehensive measurements were made of acceler- 
ated flight. 

By 1925 the first of the several reports on the full-scale deter- 
mination of the derivatives for a Bristol Fighter, BF2b,** was 


* The second edition, 1939, of this text is widely used as a basic reference. 
The dimensional notation is used. Bairstow’s original work in stability 
was so prominent that the basic approach is often referred to as that of 
“Bryan and Bairstow.” 


published by the British. These tests required several years to 
complete and are unique in that at no time since have such com- 
plete measurements of the lateral derivatives been accomplished 
full scale. With the exception of the determination of the lon- 
gitudinal derivatives for a light airplane at the N.A.C.A.* and 
isolated tests in Germany (Dornier) and elsewhere, no further 
flight tests for the complete dynamic derivatives were made. AI- 
though the agreement in unstalled flight between theory, wind 
tunnel, and full scale was at that time (14 years ago) reasonable— 
particularly for power-off conditions—the extended use of classi- 
cal dynamics in research and design has been based on little more 
full-scale flight proof than that noted. 

A dynamic problem of great complexity is that of the stall. 
In 1920, coincidental with research in the unstalled régime, the 
British launched a determined attack on this problem, making use 
of all wind-tunnel, theoretical, and flight resources (with results 
that were even less known and used by the practicing engineer 
than those from the unstalled investigations). In the course of 
this work, extensive calculations were performed in the step-by- 
step integration of the equations, thus enabling the motion to be 
followed for arbitrary control displacements. 

Tests*~* were performed in stalled glides at angles of attack 
of up to 50° which indicated the superiority of the large vertical 
tail at the stall, the effectiveness of slots and interceptors as con- 
trol, etc. While this overall program‘ gave a basic understand- 
ing of stall dynamics, much remains to be learned. Although the 
British assumed the existence of the conventional velocity deriva- 
tives, it now seems apparent that the classical approach is inade- 
quate for completely explaining even the initial divergence. In 
this country, with the exception of Koppen’s investigations, the 
problem of control at the stall has largely been thought of in 
terms of basic wing flow breakdown. The problem is a dynamic 
one, however, involving the a‘rcraft as a whole and is becoming 
more important with extreme plan form and airfoil shapes, which 
for performance reasons can be subject to but limited modifica- 
tion. 


Research in the Decade Prior to the War 


Since about 1930 the course of airplane dynamics in this 
country has been widely and continuously influenced by the re- 
search and teachings of Otto C. Koppen.*! His theoretical in- 
vestigations, wind-tunnel work (oscillators), and achievements 
in airplane design are now well known. For years his course in 
stability and control at M.I.T. was unique in its treatment of the 
complete dynamics, and many current trends in design and re- 
search may be directly traced to his work as an educator. The 
elevator forcing function technique, with which the bulk of the 
work of the present paper is concerned, had earlier been suggested 
by Koppen, and it was through his influence that the investiga- 
tions were initially directed at the longitudinal motion. 

In 1935, Mueller‘? devised one of the first simulators, an elec- 
trical device for solving the classical longitudinal stabiiity equa- 
tions. Unfortunately, a similar device for the lateral motion was 
never constructed, nor was the effect of controls introduced into 
the longitudinal machine. Probably the first application of the 
electrical engineer’s rotating vector diagram to the solution of 
airplane stability problems was presented by Mueller‘ in 1937. 
This approach deserves great attention, allowing, as it does, a 
strobotactic view of the elements of a vibration problem and their 
phase relationships. In a section of the present report, a similar 
“electrical’? method has been applied to the determination of ef- 
fective damping and spring constants from experimental flight 
data, an application suggested by I. G. Ross. In a later publica- 
tion, Mueller summarized the determination of dynamic data 
from scale-model airplanes. 

An extended series of wind-tunnel tests was initiated by the 
N.A.C.A. in 1932 for the evaluation of lateral control devices.*! 
Subsequent research indicated that the quasi-static rolling criteria 
used were inadequate, and a method based upon a step-by-step 
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solution of the classical equations was devised by Jones and 
Weick,* 1936, and applied in a final evaluation“ in 1937. Paral- 
leling the wind-tunnel research, an intensive flight program, utiliz- 
ing various devices and aimed principally at control effectiveness, 
was initiated on an F-22 airplane. These remarkable flight tests 
under Weick, Gough, Soule, McAvoy,* and others, led to first- 
hand knowledge of spoiler lag, two-control operation, effect of 
excessive dihedral, etc. The experimental methods themselves 
were widely copied by the industry. Extensive tests of the F-22 
in the full-scale tunnel were made in conjunction with the flight 
program. The backlog of experience from these tests appears to 
have been influential in Gilruth’s later flying qualities develop- 
ments. 

Zimmerman,** in two widely known N.A.C.A. reports, pre- 
sented a complete engineering analysis of the power-off longitu- 
dinal and lateral stability, with charts for design use. The ap- 
proach was entirely classical and followed in part the earlier work 
of Gates and Koppen but with an original and comprehensive ex- 
amination of the phugoid. They brought to the attention of the 
American engineer the great advantage of the nondimensional 
notation. In addition to Zimmerman’s work, the N.A.C.A. 
continued its theoretical and wind-tunnel researches into the 
evaluation of the individual derivatives, and by 1942 a large 
mass of information was available to the designer for the estima- 
tion of these constants under normal flight conditions. Refer- 
ence*® °° js here made to but a typical part of this expanding liter- 
ature. Many specific applications, such as to two-control opera- 
tion, were also considered. 

As early as 1930 Bryant and Williams*! had suggested the use 
of the method of operators for calculations of the disturbed mo- 
tion of an airplane. A simplified, and now widely used, applica- 
tion of the method was given by Jones*? in 1936, who later fol- 
lowed with publications on allied research. The method has 
since seen considerable use in the design of unconventional air- 
craft (such as flying wings), in the computation of maneuver loads 
on tail surfaces by Perkins®* and others, and in automatic control 
applications. An improvement in the graphical solution of the 
Duhamel integral has recently been suggested by Laitone* and is 
specifically referred to in the body of this report. The intro- 
duction of the method of operators to dynamic stability and con- 
trol calculations represents one of the few advances made in the 
mathematical treatment of the subject since the time of Bryan. 

With the exception of work referred to elsewhere, European 
progress in dynamic stability and control in the 10 years preced- 
ing the war was not great. In 1930, as the result of directional 
hunting on Schneider Cup racers, the British investigated® the 
directional stability of high-speed aircraft from a purely classical 
approach, a work significant in the fact that similar difficulties 
have since been encountered on dive-bombing and fighter air- 
craft and have since been more thoroughly investigated. The ef- 
fect on stability of high-wing loadings had been estimated through 
the medium of the nondimensional equations, and attention was 
being directed to the behavior of unconventional types (tailless, 
autogyro, etc.). In 1934, Jones’s? review of progress was presented 
in Durand’s Aerodynamic Theory, which is still the most widely 
used reference text. In the late thirties, the trend of British re- 
search, subjected to American influence, swung far toward appli- 
cation, with emphasis on statics. German research had, of 
course, accelerated at an unprecedented pace. 

A powerful experimental method of studying dynamic be- 
havior came into common use in Great Britain in the early thir- 
ties. This was the testing of dynamically similar models, which 
is basically older than the airplane itself (Langley’s !/,-scale aero- 
drome). Direct observation of these models gave information 
that could, through experience, be correlated with full-scale be- 
havior. In recent years the method has been applied to the hy- 
drodynamic characteristics of seaplanes by Coombes, Perring, and 
Johnston® and introduced into this country by Stout, Shoe- 
maker, Truscott, and others. The preliminary tests in the 


N.A.C.A.’s first free-spinning tunnel were presented by Zimmer- 
man*” (1936), who later originated the free-flight tunnel. At 
about this same time the N.A.C.A. gust tunnel for studying dy- 
namic loads was placed in operation. While the hydrodynamic 
porpoising and the free-spinning model tests have substantially 
replaced, in their respective fields, other methods of design, it is 
unlikely that the free-flight and gust tunnels will ever reach such 
a degree of sufficiency, dealing as they do with problems through- 
out the flight range. In its present form the free-flight tunnel 
has demonstrated its remarkable capabilities in coping with a 
wide range of dynamic problems®*® (on coaventional, unconven- 
tional, powered, stable, uistable, and towed models), the refer- 
ence cited being but typical. Though checks have been made with 
the full scale, its limitations remain as those of scale, smal! devia- 
tions from straight flight and the physical size of the model which 
pose serious model design problems. In addition to the free- 
flight model, many proposals have been put forth for semicon- 
strained models, such as a model attached only to a drag wire. 
No outstanding results with this type of constraint are known. 

One step removed from the free-flight tunnel is the free-flying 
radio, or otherwise remotely coatrolled, model exploited by Zim- 
mermazn and, more recently, exhaustively developed by Stout," 
which seems emimently suited to three-dimensional maneuvering 
problems. A further development of the flying mode! is the use 
of small man-carrying aircraft to simulate larger, heavier, or 
higher performance machines, a scheme that has been used with 
considerable success by Short Brothers (Great Britain), Martin, 
Boeing, Chance Vought, and others, with wide variations in the 
degree of similitude. This technique of flying man-carrying 
models has come under discussion relative to missile dynamics. 

In the early thirties, attention was directed to the problem of 
the desirable flight characteristics of the manually controlled air- 
plane, a difficult phase of the picture since it involved a consoli- 
dation of practical pilot experience with engineering fact. This 
had never previously received adequate attention. In about 
1935, E. T. Allen observed that large German aircraft possessed 
superior handling qualities, a trend that seemed hardly fortuitous. 
Soule*! (1936) compared the longitudinal dynamic stability of 
several airplanes with pilot observations, which finally eliminated 
the phugoid from consideration as an important flying quality, 
and about the same time Warner set down the handling specifica- 
tion for the original DC-4. By 1940, Soule had reported on his 
preliminary flying qualities investigation (an outgrowth of the 
DC-4 specification), and the work was well under way with Gil- 
ruth’s and Turner’s work of the same year. Gilruth’s®? Require- 
ments for Satisfactory Flying Qualities of Airplanes (1941) became 
a basic reference on the subject. Following this work, flying 
qualities requirements have been promulgated for military air- 
planes through the efforts of Perkins,®* ** Abzug,® and others, 
and complete tests have been performed on over 30) American 
and foreign military aircraft of all types. This impressive array 
of full-scale information has gone far in orienting effort in the 
stability and control field. Further ramifications have been 
studied by Phillips,** Kaytan, etc., and the entire subject (with 
test methods) was reviewed with the industry and services in a 
series of conferences in October, 1944, at N.A.C.A. Similar 
specifications are now available from British and German experi- 
ence. In all of these specifications, a number of dynamic and 
quasi-static requirements are included which may be modified or 
appended as new experience aad methods become available. 
The serious development of specifications like these for civil air- 
craft use may result in stringent requirements relative to rough- 
air behavior and performance on automatic pilot. The entire 
subject of desirable flying characteristics is, of course, directly 
related to research in dynamics whose function is to further the 
attainment of information that will enable the design of aircraft 
for predetermined characteristics. 

A flight method for the numerical determination of the slope 
of the static pitching moment, though not required by the flying 
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qualities specifications, has been developed to a point of practical 
utility for comparing wind tunnel with full scale, and in specifying 
safe operating c.g. limits. This method, most widely used by 
Boeing,” yields total moment change with speed, hence having 
dynamic implications. Known as the ‘effective weight moment 
method,’ it is not new in principle. 

A new tool for the measurement of the derivatives, and the 
study of dynamic phenomena in general, is the N.A.C.A.’s 
stability and control tunnel. By placing the model in a curved or 
rotating air stream, many of the disadvantages of previous dy- 
namic balances are avoided, and it should be possible to conduct 
tests at higher Reynolds Numbers. It is understood that no 
dynamic data are yet available from this tunnel. 

By 1940 a modern and comprehensive analysis of free-control 
stability was completed by Jones and Cohen.® 


Recent American and British Trends 


During the war and up to the present time, the major effort in 
dynamic research has been in the direction of application. This 
effort, largely based on established methods, has served to reduce 
partially the reluctance of designers to make use of the fund of 
dynamic information already available. (The experience of A. C. 
Reed and the author, in connection with the design of a prone 
position flying wing fighter, is illustrative of the necessity for re- 
sorting to fundamental calculations of the motion, since the mag- 
nitudes of the static stabilities were not only different from 
conventional designs but could not be judged in terms of conven- 
tional design experience. Others engaged in the design of un- 
conventional aircraft have faced like problems and have made 
computations of the motion resulting from control displacement 
by step-by-step integrations, Duhamel’s integral, or straight 
solution.) Among the many practical problems to which re- 
search has been directed during the war years are the damping of 
the lateral oscillation, aileron roll and stick force characteristics, 
the effect of friction in controls, wing-fuselage interference 
effect on the derivatives, and power effects, to cite but a few. 
The large amount of applied dynamic research still in prog- 
ress in the large research laboratories defies ready classifica- 
tion. 

In addition to applied research, some attention has been di- 
rected to more fundamental dynamic considerations. A portion 
of this work relates to the determination of the derivatives (for 
small disturbances) for various types, and degrees of rapidity, of 
motion. It will be appreciated that for extremely slow motion 
the simple Bryan derivatives depending only upon the present 
value of the velocities are probably satisfactory; on the other 
hand, for a motion of the rapidity of flutter, the aerodynamic 
forces do not follow such simple laws, and resort must be had to 
theoretically calculated flutter derivatives. It is now believed 
that even for motion in the vicinity of the longitudinal short 
period, the, simpler form of derivative may be inadequate. This 
entire subject, linked as it is with flutter research, represents a 
broad field for investigation for many years and is mentioned in 
connection with the investigations of the present report. In re- 
lation to the domain of dynamic stability and control, reference 
is made to the work of Temple, Jones, Duncan, Cicala, Reid, and 
Vincenti.“-78 Furthermore, it is becoming clear that for certain 
wing configurations, such as high sweep, it is no longer possible to 
estimate even the simplest type of derivative from present aero- 
dynamic theory. Effort in both Great Britain’ and this country 
is being directed toward further span load analysis. 

Since about 1940, Mitchell, Gandy, Bryant,%-7 and others in 
Great Britain have labored upon the calculation of airplane 
response characteristics for wide ranges of design parameters. 
The case for these extensive calculations has been well stated by 
Bryant,”4 Gates,”* and Mitchell and Thorpe.”? The response cal- 
culations have been made for various types of applied disturb- 
ances, usually the simple step. The calculations were performed 
on the differential analyzer, which represents the first time this 


tool has been put to such extensive use in the interests of airplane 
dynamic stability and control. 

Application of the differential analyzer is representative of a 
general trend toward calculators and simulators of all sorts, 
The situation regarding these various devices in stability and 
control research has recently been summarized by Blaschke.* 
Since in all cases the motion is expressible in some type of differ- 
ential equation, the advantage attendant to the use of such de- 
vices in stability and control is obvious. The reference alluded 
to® not only examines the principles of available machines but 
contains a complete bibliography of the literature. The use of 
such equipment should go far in removing the stigma that has so 
often dogged stability calculations in the past: however, the 
main problem remains that of providing the basic aerodynamic 
data. The devices now under examination vary from simple 
computers to elaborate ‘“‘Link Trainer’’ types into which arbi- 
trary control motions may be continually fed. In addition to 
machines, research has been directed to the solution of the dif- 
ferential equations by tabular, matrix, and graphical means. 

An electrical balance has been proposed by Bratt and Wright,*! 
the measuring elements being of the Eastman type (University 
cf Washington wind tunnel) in which the support reactions 
would be measured during forced oscillation of the model. The 
balance should have utility in both stability and flutter work. 
(The author, in 1942, proposed a strain-gage type of dynamic 
balance for the high-speed wind-tunnel measurement of deriva- 
tives on small models. The realistic design of dynamic balances 
for high-speed work should, it is now believed, follow further dy- 
namic flight tests.) The amount of dynamic research performed 
in high-speed tunnels appears to be extremely small. 

In the more generalized motion of the airplane in which the 
longitudinal and/or lateral motions are of finite magnitude, some 
progress has been made. This has been summarized by Bro- 
detsky®? in a contribution extending the earlier work of Lanches- 
ter and Hopf. Of course, the finite motions of airplanes have 
been considered in specific problems (landing flare, dive pullout, 
stall, etc.) for many years, but in these cases resort was invariably 
made to step-by-step integrations or drastic simplifying assump- 
tions. 

Some results have been obtained in both England and this 
country in the measurement of atmospheric turbulence, which is 
a necessary preliminary to the dynamic design of airplanes with 
improved rough-air flying characteristics. Progress has been 
made in the use of airplanes as gust meters and will be enhanced 
by better methods of specifying the basic smooth-air characteris- 
tics. In the author’s opinion, the need for rough-air response 
measurements and atmospheric turbulence studies is well sum- 
marized by Gates and Perring.§* This is a large and important 
field for research. 


Automatic Control Research 


It is not within the province of this report to examine the back- 
ground of design and analysis used in the automatic control field 
other than to emphasize again that a lack of engineering data on 
airplane dynamics has been a limiting factor in rational design. 
In general, it appears that the many automatic control studies 
that rest on the classical airplane stability equations have been 
of little use in design; conversely, most proposed simple analyses 
have failed to approximate the aerodynamic facts. In this di- 
lemma the automatic control designer has proceeded largely on an 
empirical course with a minimum of professional aerodynamic as- 
sistance. With his electrical background, he has frequently had, 
however, a better understanding of the dynamic characteristics 
of a particular airplane than its own design staff. 


Resume of .German Advances from Discussion with 
German Scientists 


From discussion with German scientists formally engaged in 
stability and control research for airplanes and missiles, it appears 
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that important advances have been made in the general handling 
of aircraft dynamics, particularly in connection with automatic 
control where refined methods were in general use. It was the 
practice in the larger aircraft firms to perform complete dynamic 
analyses for each new aircraft. 

The basic approach to dynamic problems was essentially classi- 
cal. As a result of a comparison of the different reference axes by 
Rautenberg, it was concluded that for many purposes axes fixed 
in the flight path were most desirable. The limitations of the 
derivatives at higher frequencies were recognized. The restric- 
tions in wind-tunnel measurements of the derivatives were well 
recognized and other methods of research were sought, one of 
With this de- 
vice, sizable ‘‘aerodynamic forces’’ could be obtained in compari- 
son to the centrifugal forces involved which should improve the 
measurement accuracy over the conventional (air) whirling arm. 
A complete dynamic strain gage balance was designed (Thiel), 
Lut it was not built when calculations of its basic accuracy led to 
the belief that the tares would be excessive. Sufficient wind- 
tuanel measurements of N, with running propeller were made to 
eusure the accurate calculation of this derivative. Research di- 
rected toward the experimental measurement of NV, on the rolling 
balance has been reported by Thiel. 

As in other countries, only a small amount of full-scale airplane 
research was performed (other than in connection with automatic 
control design). Some of this work has been reported by Schmidt 
and Stiess, was done on a Dornier DO-117, and was largely con- 
fined to the static lateral derivatives Y,, N,, L,. Other flight 
measurements have been reported by Wenk and Doetsch of 
D.V.L. Insofar as the author could determine, no sinusoidal 
control forcing function tests were made in flight, the closest ap- 
proach to this technique being that used by von Scalippe (Junk- 
ers, Dessau) in flight flutter tests. 

In the field of automatic control, many unique possibilities were 
investigated. Among these were the possibilities of automatic 
control in counteracting the instabilities found at high Mach 
Numbers, the improvement of rough-air flying characteristics, 
and the increase in damping about a particular axis of an air- 
plane otherwise flown manually. Some typical references to the 
automatic. control dynamic work are given,*4~® it being appre- 
ciated that these are chosen almost at random from several hun- 
dred allied reports. 

Sufficient experimental work was done with the stall to indi- 
cate its primary dependence on angle of yaw and rolling velocity 
and to determine that these would be controlling variables 
for automatic control of the stall. 

Control step functions, or ‘‘personal gusts,’” were used in mis- 
sile flight test, with many variations such as the sudden deflection 
of two controls to give the missile a spiral motion from which 
various information could be computed. Some of the advances 
made in the dynamics of torpedoes -appear to be transferable to 
aerodynamic cases. 

All in all, there seems to be little question that the Germans 
fully appreciated the potentialities of dynamics and the utter 
dependence of future airplane and missile design upon it. 
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Letter to 


Dear Sir: 

In the excellent summary of the results from the spark photog- 
raphy range in their relation to aerodynamics (Some Ballistic 
Contributions to Aerodynamics, A. C. Charters, JOURNAL OF THE 
AERONAUTICAL SCIENCES, Vol. 14, p. 155, March, 1947), Dr. 
Charters makes reference to some work* by the present writer 
on head-base interlocking effects on the drag coefficient. The 
statement is made that this analysis shows the drag to be roughly 
independent of body length for bodies one to three calibers long. 
I feel obliged to point out, however, that although the analysis 
referred to was characterized by a very limited amount of data, 
the results indicated a systematic dependence of drag coefficient 
upon body length. That which makes the results worthy of 
note was the direction of the effect, the projectile with the longer 
body exhibiting the smaller value of the drag coefficient. 

The experimental data consist of drag coefficient measures on 
19 projectiles with body length ranging from 1.84 to 2.00 in., 
and on four projectiles with body length 2.79 in., using the ex- 
perimental values of Ap for the 19 ‘‘short’”’ projectiles and 
smoothing via the ‘‘Q’’ factor described in Charters’ paper, Kp 
values for the Mach Numbers at which the four ‘‘long”’ projec- 
tiles were tested and computed. These were compared with the 
observed results and corrected for yaw on the ‘‘long”’ projectiles. 


* Thomas, R, N., Ballistic Research Laboratory Report 543. 


the Editor 


M Kp (Short) Kp (Long) 
1.848 0.146 0.143 
1.899 0.144 0.141 
2.434 0.121 0.120 
2.527 0.118 0.115 


Obviously, the material is hardly complete enough to warrant 
quantitative conclusion. Nonetheless, in view of the internal 
accuracy of the results (about 1 per cent), the effect must be con- 
sidered real. The models were identical save for body length. 
The longer bodies cannot fail to have an added contribution to 
Kp arising from the added frictional drag. (Kp is based on the 
cross-sectional area.) Thus, the result must be taken to indicate 
some effect of added body length upon conditions at the base. 


A tentative suggestion, but only tentative, arises from reflec- 
tion upon the properties of a gas expanding around a corner. 
For the same free-stream pressure, a decrease in the free-stream 
velocity, holding fast the angle of expansion at the corner, corre- 
sponds to an increase in pressure in the stream after the expan- 
sion. Thus, for a decrease in the velocity of the gas along the 
projectile body, accompanying in any way the added body length, 
we might expect to presage a decrease in drag. This is the effect 
observed. 


RICHARD N. THOMAS 
Harvard College Observatory 
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Motion of an Unarticulated Helicopter Blade’ 


JONATHAN WINSONt 
California Institute of Technology 


SUMMARY 


This paper investigates the bending oscillation perpendicular 


. to the plane of rotation of the unarticulated helicopter blade. 


The motion is comparable to the flapping of the hinged rotor 


blade. 
Two differential equations of motion are written in terms of 


two of the natural modes of the blades. The solutions of the 
equations are a pair of infinite Fourier series in terms of the posi- 
tion of the blade in azimuth. The series are approximated by 
their constant and first harmonic terms. Six simultaneous linear 
algebraic equations written in terms of flight and blade param- 
eters give the six Fourier coefficients. A sample problem is car- 
ried through, beginning with the physical characteristics of a 
blade and yielding the bending motion at a specific flight condi- 


tion. 
The method can be expanded both to include a larger number 


of natural modes and to achieve a closer approximation to the 
infinite series solutions. There are included a sufficient num- 
ber of variables to describe blades of current design. 


INTRODUCTION 


6 bes FLUCTUATING AERODYNAMIC FORCES that cause 
flapping in the articulated rotor bend the blades of 
the unarticulated or rigid rotor in a vertical plane. 
This pure bending oscillation is treated here. Coupled 
motions are not considered. They are effectively elimi- 
nated by the assumption that the mass center, aero- 
dynamic center, and elastic axis of the blade are at a 
point at each cross section and that the elastic axis is 
a straight line. 

The method that is presented predicts the motion of 
the blade under the condition of steady forward flight 
of the aircraft. It is the direct counterpart of the solu- 
tion for the flapping motion of the articulated rotor and 
is the first step in the rational analysis of the rigid rotor, 
since the flapping equations are the necessary beginning 
of hinged rotor analysis. 

The immediate external effect of the bending of the 
blade is to cause fluctuating moments on the rotor hub. 
These moments are the most obvious cause of rigid 
rotor vibration and a basis for the investigation of such 
vibration is one of the most important functions of the 
present analysis. In addition, knowledge of the bend- 
ing motion yields such information as the fatigue cycle 
of the blades, the lag angle required in the aircraft con- 
trol system, and the maximum deflection of the 


blades. 
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OUTLINE 


The work is divided into the following five sections: 


I. The Derivation of the Bending Forces That Exist 
in Forward Flight Along a Blade. 

The differential force is expressed as a three- 
harmonic Fourier series in terms of the azimuth posi- 
tion of the blade. The coefficient of each term is a 
power series in blade radius. The aerodynamic 
theory is based on work by Glauert.! 


II. The Calculation of the Natural Modes and Fre- 
quencies of the Blade. 

(a) The method of Myklestad* is summarized. 

By considering the blade made up of a finite number 

of discrete masses, a tabular method is developed 


which allows the calculation of the natural modes and | 


frequencies of the blade in its centrifugal field. 
(b) A sample blade is considered and the three 
lowest modes and frequencies are found. 


Ill. The Derivation and Solutions of the Equations 
of Motion. 

LaGrange’s equation is applied to the system using 
as generalized coordinates arbitrary functions of time 
which multiply the two lowest modes. The gener- 
alized force is derived in terms of the aerodynamic 
force. There result two simultaneous second order 
linear differential equations with variable coefficients. 
The solution of the equations is a pair of infinite 
Fourier series in terms of the azimuth angle. The 
series are approximated by their constant and first 
harmonic terms. The final solution for the blade 
motion is given by a set of six simultaneous algebraic 
equations in terms of the blade and flight parameters. 
This app‘ication of the LaGrangian was suggested by 
work of M. A. Biot.* 


IV. A Numerical Example. 
The calculation of the motion of the sample blade 
_is carried out and plotted at a condition of forward 


flight at 100 m.p.h. 


V. Discussion. 
A discussion of the assumptions involved and com- 


ment on the application of the method are given. 


SECTION I.—THE AERODYNAMIC FORCES ON A 
ROTATING BLADE 


Nomenclature 
R_ = radius of the blade : 
r = coordinate distance along the blade span 


fo 


shaft p 


91, 2 


Thrust 


Fig. 
element 
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MOTION OF AN UNARTICULATED HELICOPTER BLADE 


y = coordinate distance perpendicular to blade span, defin- 
ing position of deflected blade with respect to its un- 
bent position 

c = chord length at any blade section 

& = chord length at the root of any tapered portion of the 

blade 

coefficient describing the rate of decrease of chord with 

span for any tapered portion of the blade. For a 

blade linearly tapered from root to tip: c = @ — 

re(r/R) 
The plane of the rotor disc is the plane perpendicular to rotor 
shaft passed through the rotor hub. 

% = geometrical blade setting of the blade at its root meas- 
ured from the line of no lift of the airfoil section 

Ke = coefficient describing the rate of change of blade angle 
with span in a linearly twisted blade. Thus, @ = 
6o — Ke(r/R) indicates linear washout from root to 


Ke 


tip 

6 = geometrical blade angle at any station along the span, 
in general a function of r and time 

y = azimuth position of the blade measured in the direc- 


tion of rotor rotation. The blade in its downwind 
position is at zero angle of y 


u = velocity of flow induced through the rotor disc, as- 
sumed constant over the disc 

V; = component of the resultant velocity of the air relative 
to the blade in a plane perpendicular to the blade 
span 


Vj, = component of V; perpendicular to the blade span and 
lying in the plane formed by the blade span and the 
rotor axis of rotation. 
V;z = component of V; which is perpendicular to both V;, 
and the blade span 
dT = thrust component of resultant aerodynamic force on 
segment of blade dr. dT is parallel to V;y 
6, 62 = constants that are used to define the cyclic control im- 
posed on the rotor. Arbitrary feathering control 
produced by a swashplate may be expressed as 
Ad = 6, sin + cos 

= absolute angle of attack of the blade element 


Qa 

¢ = the induced angle; the angle between V, and V;z 

w = speed of angular rotation of the rotor P 

a = angle of attack of the rotor disc; the angle between the 
flight path and the rotor disc, positive when the disc 
is tilted forward 

p = density of the atmosphere 

a. = slope of the lift curve for infinite aspect ratio 

C, = coefficient of lift of the blade element 

t = time 

Thrust Loading of the Blade 


Fig. 1 illustrates the velocities and forces at a blade 
element. The following assumptions are made: 


Eq. (2) through Eq. (9) are substituted in Eq. (1). 


Dw 


Fic. 1. Velocities and forces at a blade element. 


(1) The angle ¢ is small so that ¢ = V,,/V,, and 
cos ¢ = l. 

(2) lV, = |Vrel. 

(3) The blade sections operate below the stall. 

(4) «is constant over the disc. 

(5) The component of blade element drag in the 
direction of thrust is negligible so that d(Lift) = dT. 

(6) Oy/dr is small so that sin [arctan (Oy/Or)] = 


oy/Or. 
Then: 
dT = (p/2)C,V,°c dr (1) 
CL = (2) 
a= 6 (V,y/ Vez) (3) 
= — Ker + sin + (4) 
c = @ — K.(r/R) (5) 
= wr + Veosisiny (6) 
Vy, = Vsni+ut (7) 
or: 

Let there be defined: 
u = Vcosi/wR (8) 
vy = (Vsini + u)/wR (9) 


There result three parts to the expression for d7: one inde- 


pendent of blade shape, one dependent on blade shape, and one dependent on the rate of change of the shape with 


time. 


series in. The final expression for dT is: 


Each part is expressed as a Fourier seriesin y. The coefficients of the Fourier series are arranged as power 


dT = dr[(Ay + Ai sin + Az cos + Az sin 2y + A, cos + As sin 3y + Ae cos 3y) + 


(Oy/Or) (Bz cos + Bs sin 2p) + (1/w)(Oy/Ot)(Co + Ci sin 


where 


(10) 
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= |+ + OK.) + coR(Om — »| + 
(p/2)a, 
K Kop? ook, K.K 
(«ut + K + Kw) + r(- + (11) 
=I CoO mas Conv R? ri— —K + Ao + K wvR + 
(p/2)a.. 4 4 
2K.K. 
— 2coKou — 2K + — (12) 
R R 
(p/2)a., (; + (- q + r (<a) + r(- (13) 
A;/(p/2)ae = + K (14) 
A Co KC K Kop? 
= ( +r ) + r(- + K dun) (15) 
As/(p/2)a,, = + (16) 
Ae/(9/2)a,, = + r(1/sK Oou?R) (17) 
Bo/(p/2)a,, = r(—coR) + (18) 
B;/(p/2)a,, = (— cou?R?/2) + r(K u?R/2) (19) 
Cé/(p/2)a,, = r(— co) + r°(K,/R) (20) 
Ci/(p/2)a., = (— couR) + r(K (21) 


Section ITA.—TueE NaTuRAL MODES AND 
FREQUENCIES OF THE ROTATING BLADE 


The Myklestad method for the calculation of the 
natural modes and frequencies of a rotating blade may 
be found in reference 2. A word summary of the 
method is given here. j 

The theory of normal coordinates indicates that, at a 
natural frequency, the rotating cantilever beam consid- 
ered will vibrate in such a manner that the following 
conditions are satisfied : 

(1) The transverse motion of all points on the beam 
will be harmonic and in phase. 

(2) The motion will exist without application of 
external forcing functions. The external forces acting 
on a point mass will be only those of inertia and centrif- 
ugal force. 

(3) The motion will be consistent with the elastic 
characteristics of the beam. 

(4) The boundary conditions of the system will be 
satisfied. They are: (a) deflection at the root = 0 and 
(b) slope at the root = 0. 

By a proper mathematical statement of the four 
conditions, it is possible to find both the frequencies at 
which this motion will exist and the mode of motion at 
each frequency. 

More specifically, the system is assumed to consist of 
a finite number of discrete masses connected by a mass- 
less elastic beam. The four points mentioned above 
are imposed with the exception of boundary condition 


4b. This boundary condition is made a function of an 
assumed natural frequency of the blade. When the 
assumed natural frequency is the true natural fre- 
quency, the boundary condition is satisfied; that is, the 
slope at the root of the blade is zero. Thus, a plot of 
assumed natural frequency against calculated value of 
the slope of the blade at its root yields as many of the 
natural frequencies as are desired. The mode associ- 
ated with each frequency is found with little additional 
calculation. 


SEcTION IIB.—NATURAL FREQUENCIES AND MODES 
OF THE SAMPLE BLADE 


The sample blade considered is shown in Fig. 2. Its 
main structural member is a step-tapered hollow steel 
spar. So little stiffness in vertical bending is added by 


ACTUAL BLADE 


FINITE _MASS_APPROX/MATION 


4 45 7™ 102 129 181 202 WCHES 
Fic. 2. The sample blade. 
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TABLE 1 
Physical Characteristics of the Sample Blade 
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Station, In. Spar Size, In. Weight, Lbs. 
O0- 12 3.50 X 0.263 9.09 
12- 30 3.25 xX 0.224 10.86 
30- 45 3.00 X 0.200 9.08 
45-— 60 2.75 X 0.183 9.21 
60- 74 2.50 xX 0.170 8.02 
74— 90 2.25 X 0.148 7.80 
90-102 2.00 X 0.131 5.02 
102-114 1.75 0.117 4.30 
114-126 1.50 X 0.109 3.79 
126-145 1.25 X 0.104 5.17 
145-169 1.00 X 0.104 5.66 
169-193 0.875 XK 0.102 4.90 
193-210 0.750 X 0.102 3.05 

85.95 Ibs. 

TABLE 2 


Finite Mass Approximation of the Blade 


Station, Concentrated Mass, 
In. Lbs.Sec.?/In. 
14 0.0517 
45 0.0474 
74 0.0410 
102 0.0241 
129 0.0232 
157 0.0147 
181 0.0127 
202 0.00790 


the covering that the spar is considered to be the only 
elastic member present. The weights of the covering 
and ribs are taken into account. 

The essential characteristics of the blade are given 
in Table 1. Spar size is written as diameter X wall 
thickness. Table 2 lists the system of masses chosen 
to represent the blade. Fig. 3 is the plot of assumed 
frequency versus slope at the root of the blade which 
was used to find the three lowest natural frequencies. 
Fig. 4 is a plot of the modes associated with these fre- 
quencies. 


SEcTION III.—MOTION OF THE BLADE 


The motion of the blade under the applied aerody- 
namic forces will now be found. In this section the sys- 
tem is considered continuous rather than made up of a 
finite number of discrete masses. 

The motion, in general, will consist of some combina- 
tion of the natural modes. Only the effects of the low- 
est two natural modes are considered here. Section V 
discusses this point. 

The normal modes (normalized by taking unit deflec- 
tion at the blade tip) may be described by power series: 


y = (22 
and 
y = y2(r) (23) 
The complete motion may be expressed as: 
Y= + (24) 


where 


= gilt) (25) 


O10; 


SLOPE AT ROOT - KADIANS 


NATURAL FREQUENCY - RAD SEC. 


Fic. 3. Assumed natural frequency versus slope at the root of 
the blade. 
of 364 
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Fic. 4. The three lowest natural modes. , 


gz = go(t) (26) 


gi and q, are taken as the generalized coordinates in 
an application of the LaGrange equation. The gen- 
eralized forces are found from the expressions for the 
aerodynamic force. 


LaGrange’s Equation: 
(d/dT)(0T/0g;) + (0U/0q:) = Q: (27) 


where 
T = the kinetic energy of the system 
U = the potential energy of the system 
Q,; = the generalized force, as defined later 
gi = dq;/dt 


The Kinetic and Potential Energies: 
T = 1/, dm (28) 


where dm is the differential mass. 
(24) in Eq. (28) and squaring, 


T = dm + '/2 dm + 
Because of the orthogonality of the principal modes, 


dm = 0 (30) 


Substituting Eq. 


| 
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Thus, 


The theory of normal coordinates also indicates that 


the potential energy is simply related to the kinetic © 


energy through the natural frequencies and may be 
written as follows: 


where 2; and 2) are the lowest and next highest natural 
frequencies of the blade. 


The Generalized Force: 


Q, is defined as 61V/6q;, where 6W is the work done 
by the external forces on the system during a small 
change in the coordinate g; The small change in co- 
ordinate is indicated by 6g;. In this case the only ex- 
ternal load is that of the air. Its magnitude is dT in 
length dr of the span. 


bW = JS dT by (33) 
over the blade, where éy occurs because of 6q;. 
by = (34) 


Substituting Eq. (34) in Eq. (33), 
bW = dT = 6g: dT y, ~— (35) 


Thus, 
Q,= S,* aT y, (36) 


The equation of motion for the coordinate q; is now 
written in the following manner: 

(1) Eq. (24) is substituted in Eq. (10). 

(2) The resulting expression for dT is substituted 
in Eq. (36), 7 being set equal to 1. 

(3) Eqs. (31) and (32) and the expression for Q; 
found in step 2 above are substituted in Eq. (27), the 
LaGrange equation. 

(4) The symbolism y,’ = dy,/dr is introduced and 
the following constants are defined: 


= Serv? dm 

As = dr 
dr 
hr = y:Aadr 
= 
Xo = y As dr 
Aw = Be dr 
Au = yin’ Bs dr 
Me = So. dr 
AB = Bs dr 
Au = dr 
Ms = dr 
Aw = dr 
Av = dr 


1947 


A change of variable is made. Recognizing that y = 
wt, the final equation is written in terms of y. In terms 
of this variable the differential equation of motion for 
the coordinate g, becomes: 


ad’ d 
Ai Aus — Aas sin 
2? d 
a Aw cos duu sin 24 ) + x 


(— Aw — Ar7sin + cos — Ais sin 2y) = 
As + Ay sin + A; cos + As sin + 
dz cos + As sin + Ay cos (37) 


A similar set of substitutions, in this case with 7 being 
set equal to 2, leads to the equation for the variable q. 
The set of constants are: 


de = _ dm 
=> o dr 
ha’ = dr 
de’ = wAs dr 


7’ = yoA, dr 
es’ = dr 
dg’ = rAs dr 
dio’ = te. By dr 


An’ = Bs dr 
yoyo’ Bo dr 
A’ = yoyo’ Bs dr 
Au’ = y2yiCo dr 
Ms’ = dr 
ue’ = Wx, y2?Co dr 

= Cy dr 


The equation is: 


~ 
ll 


d? 
(— dw’ — Aw’ sin ¥) + 
2 d : 
2 cos y Ais’ sin 24) + x 
w dy 


(— Au’ — As’ sin Y) + Aw’ cos Y — 
Au’ sin = Az’ + Ay’ sin Y + As’ cos + 
de’ sin + Ax’ cos 2W + As’ sin + 
dy’ cos (38) 


Solution of the Equations: 


The two infinite Fourier series, which are the solutions 
of Eqs. (37) and (38), are approximated by their con- 


stant and first harmonic terms. Thus, 
= % + a sin + a2 cos (39) 
= bo + by sin + cos (40) 


Substitution of Eqs. (39) and (40) in Eqs. (37) 
and (38) and harmonic analysis leads to the follow- 
ing six algebraic equations for the six Fourier coeffi- 


cients: 
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Q;° Aus — A 


a2 + + — = Ag’ (42) 


= hy’ (44) 
Au 
+ 1— As 2+ 
2 
Ais 
+ (Ars + b = — Xs (45) 
, An’ 
Go + (as + + + 


Ais’ Q,? 


The solution of these six equations gives the motion 
of the blade at a specific flight condition within the 
approximations made. 


SECTION IV.—NUMERICAL EXAMPLE 


The motion of the sample blade at a given flight condi- 
tion will now be found. The following quantities must 
be evaluated: 


yn(r), yo(r); 1, Q; Ao, Ai, Ao, Bo, Bs, 


Co, Ci; A, do, Az, Na; Aro, Au, Ar, Aus, Au, Aus, As, 
Yaz, As’, Aa’, As’, Azo’, Ara’, Ano’, Aus’, Ara’, Aas’, Are’, 


Quantities 1, v2, V1", V2"! 


The following series are adequate for the representa- 
tion of the two modes: 


0.00354r? 
—0.009398r? + 0.00004584r* 


where 7 is in feet. Then, 


= 0.00708r 
‘= —0.018796r + 0.00018336r* 


Il 


Physical Constants: 


= 1.593 ft. 

K, = 1.039 ft. 

Ke =0 

R =17.5ft. 

w = 26 rad. per sec. 


p = 0.002378 slugs per cu.ft. 


=5.3 
dm = (0.258 — 0.0120r)dr 


Flight Condition: 


The following aerodynamic parameters are taken: 


nw = 0.304 

v = 0.0853 

6. = 0.2269 rad. 
6; = = 


This corresponds to flight at 100 m.p.h. of a ship 
using two two-bladed rotors of the type being con- 
sidered. The weight of the ship is approximately 3,000 
Ibs. 


Aerodynamic Quantities: 


From Eqs. (11) through (21), 

Ap = 0.03223 — 0.016197 + 0.002837r°? — 
0.00008489r5 

A, = —0.07972 + 0.027217 — 0.0009031r° 

Ay = 0 

B, = —0.05341r + 0.001991r? 

B; = —0.1421 + 0.005296r 

Cy = —0.01004r + 0.0003741r? 

C, = —0.05341 + 0.001990r 


Ni = 0.3414 
Ng = 0.7052 
As = 0 

= —0.2005 
Au = —0.03994 
Aye = —0.5364 
Aus = —0.2657 
Ais = —0.1003 
Ars = —0.04737 
= —0.00099 
d3’ = —0.08594 
ry’ = —0.07300 
As’ = 0 
Aro’ = —0.00204 
An’ = 0.00866 
Are’ = —0.3281 
Aus’ = —0.05793 
Au’ = —0.04737 
Aas’ = —0.00099 
Are” = —0.1886 
= —0.08088 


(2,/w)? = 1.9600 

(Q,/w)? = 10.388 

Substituting in Eqs. (41) through (46) and solving, 
there is obtained: 


ms 3 
for 
se] ll 
w* 2 is 
Ais 
(m (43) 
An’ 2 
dz + re — + 
w 
37) 
ing 
q2. 
Integrated Constants: j 
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Fic. 5. Deflection at the tip of the blade; flight at 100 m.p.h. 


a = +1.20 
a, = +1.00 
a, = —1.54 
bo = —0.03 
b, = —0.06 
b, = —0.01 


The motion is then described by 


y = (1.20 + 1.00 sin y — 1.54 cos W)(0.00354r?) + 
(—0.03 — 0.06 sin y — 0.01 cos y)(—0.009398r? + 
0.00004584r*) 
where ¢ is in feet. 
The motion of a point at the tip of the blade is plotted 
in Fig. 5. 


SECTION V.— DISCUSSION 
Aerodynamics 


The aerodynamic derivation is based on steady flow 
wing theory. Air forces are determined by instantane- 
ous values of velocity and angle of attack at each blade 
element. However, the following facts are known 
to be incompatible with the assumption of steady 
flow: 

(1) The blade element oscillates vertically with an 
amplitude that is large as compared to the chord. 

(2) The blade element changes its pitch periodi- 
cally. 

(3) The velocity at the blade element varies 
greatly. 

(4) The blade operates in a wake produced by itself, 
other blades of the rotor, and other rotor(s). 

Nevertheless, the Glauert steady flow theory! has 
been widely applied in helicopter work and has given 
results that have been found reasonable. Specifically, 
the calculation of the flapping motion of hinged blades 
has been carried out using this theory and has proved to 
be in satisfactory agreement with experiment. A more 
exact theory has not yet been developed. Points 1 and 
2 have been considered for the case of small oscillations 
by von Karman and Sears‘ and point 3 by Isaacs,5 
but no complete evaluation of the Glauert theory can 
be made on the basis of these. The one approximate cal- 


culation that can be made does not indicate a large 
error.* 


The Differential Equations 


The differential equations are not exact in that the 
modes assumed to make up the motion and their natu- 
ral frequencies are calculated for the undamped, con- 
servative system. This is the most compromising of 
the assumptions made and cannot be rigorously de- 
fended without an exact solution of the problem. It 
has some precedent in flutter work and in the assump- 
tions made in lesser degree-of-freedom systems (that 
a reasonable amount of damping does not change the 
undamped natural frequency), but the assumption is 
intuitive. Structural damping is neglected as unimpor- 
tant as compared to the aerodynamic damping. 


Solution of the Differential Equations 


The solutions of the differential equations are a 
pair of infinite Fourier series in y. The justification 
for approximating the series by constant and first har- 
monic terms is based on comparison with the flapping 
equation for hinged blades, a degenerate case of the 
equations derived here. The forcing function in both 
systems is characterized by strong first harmonic terms 
with higher harmonics of rapidly diminishing magni- 
tude.?- An exact solution of the flapping case* shows 
a corresponding, characteristic, strong first harmonic 
response with rather small higher harmonic motions. 
The inclusion of two modes in the present analysis al- 
lows any higher harmonic to be considered whose prox- 
imity to resonance emphasizes its importance. There is 
no theoretical barrier to carrying the series solutions 
out to higher harmonics. 


* Assuming small vertical oscillations, the effects of deviations 
1 and 2 are found to be small. The reduced frequency of the 
system as defined in the von Karman paper is of the order of 0.04. 
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In-Flight Icing of Highly Electrified Aircraft 


GUNN* 
Division of Physical Research, U.S. Weather Bureau 


ABSTRACT 


The suggestion that a highly electrified airplane will not ice up in 
normal icing conditions is examined. Tests with an artificial 
charger on an aircraft in flight and on electrified probes mounted 
on an experimental airplane have shown that the electrification 
has no measurable effect on aircraft icing. 


INTRODUCTION 


4 lpm SERIOUS HAZARD resulting from the formation of 
ice on flying aircraft has engaged the attention of 
aeronautical engineers and scientists for many years, 
and various suggestions have been made for its elimina- 
tion. It has long been known that the water particles 
in the atmosphere carry appreciable electric charges. 
A short calculation using available data! shows that the 
electrical charges on the droplets are sufficiently large so 
that one might reasonably expect the forces on water 
droplets, when exposed to actually observed electric 
fields, to affect appreciably the deposition of ice- 
forming particles on aircraft in flight. This possibility 
has been pointed out by several investigators. H. B. 
Barakan, of the Leningrad Institute of Meteorology, has 
published a paper? entitled ‘‘One of the Possible Causes 
of Airplane Icing” in which he suggests that the electri- 
fication of aircraft plays an important part in the icing 
of aircraft. More recently, a U.S. Patent No. 2,347,114 
has been issued on a ‘“‘Method and Means for Preventing 
Ice Formation on Surfaces” in which definite claims are 
made as to the influence of a large free charge carried 
by the airplane on the icing of the frontal surfaces and 
which particularly states that a suitable electrostatic 
charge on the surface of the airplane will prevent icing. 


RESEARCH 


The possible importance of any method for suppress- 
ing icing of aircraft and the fact that the suggestions 
appeared to have some quantitative foundation encour- 
aged the research team to explore the problem. Ex- 
perience and flight measurements have shown that the 
electrical state of water particles sprayed into a wind 
tunnel is in no way related to the charge on natural pre- 
cipitation droplets. It was a foregone conclusion, there- 
fore, that experiments looking toward a complete an- 
swer to this important problem and carried on in a wind 
tunnel would be subject to failure or to gross misinterpre- 
tation, because the water droplets in a wind tunnel would 
not be sufficiently aged to pick up equilibrium electrical 
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charges such as are encountered in actual ice-forming 
clouds. 

In view of these facts and because there were avail- 
able at Minneapolis an airplane and facilities that could 
be used without diverting them from their principal 
duty, it was decided to test the idea under actual flying 
conditions. By special apparatus, the B-25 assigned to 
the project could be charged artificially in flight up to 
250,000 volts and of either sign.* Two preliminary ex- 
perimental flights using this artificial charging were un- 
successful because the research team could not be sure 
of finding identical icing conditions on two successive 
runs. In order to avoid this difficulty, another simple 
experimental arrangement was adopted and led to a 
clear-cut result. In this experiment three rods */s in. 
in diameter and effective length of 4'/: in. were turned 
out of brass and carefully smoothed. These three 
identical rods were mounted on the top of the B-25 just 
behind the pilot in the Plexiglas window and about 4 in. 
in front of the aluminum skin. The three rods were 
provided with bosses and mounted symmetrically with 
respect to the axis of the ship and spaced from each 
other by 1'/. in. The three rods are shown in Fig. 1 
(the streamline body immediately behind the rods has 
no bearing on this experiment). A high voltage source 
was arranged so that any voltage from 0 to +10,000 
volts could be placed on the starboard outside rod, 
while, simultaneously, a voltage of from 0 to — 10,000 
volts could be placed on the port outside rod. The 
middle rod was grounded to the frame of the airplane 
and had no voltage on it. In general, the rods were 


operated in such a manner that positive potential of 
8,000 volts was placed on the starboard rod while 
— 8,000 volts was placed on the port rod. 


Fic. 1. Electrified probes mounted on B-25. 
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The airplane was flown through icing conditions 
twice with no voltage applied to any of the rods, and 
these blank runs demonstrated that the position of the 
rods was such that they acquired ice equally, as deter- 
mined by observation and measurement.f 

With voltages applied to the rods, the ship was flown 
to the top of the stratus cloud deck above which existed 
a slight temperature inversion. The electrodes had 
accumulated some ice during the climb, but this quickly 
melted clean in the warm air on top of the clouds. The 
top of the cloud deck was 3,300 ft. above mean sea level, 
and the temperature on top was 4°C. The high voltage 
was applied to the electrodes and the airplane dropped 
to 2,500 ft. into an icing condition that could be de- 
scribed as a moderate glaze in an ambient temperature 
of —3°C. The airplane remained in this cloud deck for 
12 min., and a substantial coating of ice on the elec- 
trodes was accumulated. It is estimated that the 
thickness of the ice layer on the leading edge of the rods 
was */,in. As far as could be observed no difference in 
the amount or distribution of ice on any of the three 
electrodes could be detected. Throughout the run the 
electric field on the top of the fuselage of the B-25 varied 
between 50 and 100 volts per cm. (a small value). 

Because it is recognized that the electrical effects of 
icing may depend critically upon meterological condi- 
tions, one other careful test was made and other casual 
runs have been made under different conditions. The 


t The author acknowledges the helpful cooperation of John 
Randolph Galt, Chief Pilot assigned to the Precipitation Static 
Project, in running these series of tests. 


second careful run, made on the succeeding day, was 
much like the first one except the rate of icing was very 
much lower and the artificial charger on the B-25 was 
operated in such a way that the airplane itself had no 
resultant free electrical charge. Again, the ice accre- 
tion was identical on all electrodes as far as could be 
determined. The icing under this condition is shown in 
Fig. 1. This photograph was taken almost immedi- 
ately upon the return of the airplane to the hangar. 


CONCLUSIONS 


Although we would not like to claim that a thorough 
investigation of the matter has been made, it is thought 
that our experiments demonstrate rather clearly that 
electrical effects are unimportant in determining the de- 
position of ice on aircraft. The electrical fields existing 
in the vicinity of the rod electrodes are approximately 
as large as any electric field likely to exist on an aircraft 
éven in dry snow conditions. Because of this fact and 
because ice invariably forms on the leading edge of air- 
craft, our experiments seem to show that the artificial 
electrical charging of an aircraft to suppress icing is of 
no practical interest. 
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The Ring Airfoil in Nonaxial Flow 


HERBERT S. RIBNER* 


National Advisory Committee for Aeronautics 


SUMMARY 


The ring airfoil with axis inclined to the stream is investigated 
by a lifting-line theory. The stream tube that threads the ring 
is found to deflect like a rigid cylinder. The lift is twice the lift 
of a flat elliptic wing that spans a diameter and has a quarter the 
area, and the downwash and ratio of lift to induced drag are the 
same for both. 


INTRODUCTION 


i be LIFTING-SURFACE THEORY of the ring airfoil in 
axisymmetric flow has been treated in references 1 
and 2. The pressure distribution is correspondingly 
axisymmetric, and there is no net lift. The ring air- 
foil with axis inclined to the stream does not appear to 
have been investigated. Such an airfoil experiences a 
net lift, and this lift is calculable by a simple lifting-line 
theory. The approximation is of the same order as in 
the lifting-line theory of the flat w.ng. The theory is 
presented herein. 


SYMBOLS 
V = stream velocity 
p = mass density of fluid 
R = radius of ring airfoil and of cylindrical wake 
¢ = chord of ring airfoil 
y = lateral distance from ring center 
6 = arccos y/R 
a = angle of attack 
a; = induced angle of attack 
é€ = angle of downwash 
= circulation 
L = lift 
C,’ = lift coefficient [L/(!/2)pV2S’] 
S’ = area of reference wing [(7/2) Rc] 
A = aspect rati> 
A’ = aspect ratio of reference wing (8R/mc) 


THEORY 


Assume, subject to later proof, that the inclined ring 
airfoil deflects as a rigid body the stream tube that 
flows through it. Let the angle of deflection (downwash 
at infinity) be 2e. For small inclination of the ring air- 
foil, this stream tube is essentially a circular cylinder. 
The downward deflection implies a distribution of axial 
vorticity along its boundary of strength: 


dl /Rd@ = 4Ve cos (1) 
where V is the stream velocity and @ the angle with the 
horizontal. (See Fig. 1.) 
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The flow velocity induced radially inward at points 
on the boundary far back is 2Ve sin 6. The radial in- 
duced velocity at the end of the cylinder that meets the 
ring airfoil is just half this. Thus, the induced angle of 
attack is 

a, = esin@ 
If the angle of attack of the axis is a, the geometric 
angle of attack of points along the ring airfoil is a sin 0. 
The effective angle of attack is therefore (a — e) sin 0. 

According to the basic assumption of lifting-line 
theory, 


= rVc(a — a,) 
where c is the section chord and the two-dimensional 
lift-curve slope is taken as 27. Therefore, 
T = rVc(a — e) sin 8 (2) 
The strength of the trailing vorticity is thus 
/Rdé@ = (xVc/R)(a — «€) cos (3) 


Now Eq. (3) shows the cosine distribution of vorticity 
characteristic of a rigidly deflected cylinder of air. 
This proves the original assumption. The identifica- 
tion of Eqs. (1) and (3) now leads to the determination 


of 
4Vecos 6 = (xVc/R)(a — cos 
a/[1 + (4R/zc)] (4) 


Returning to Eq. (2), the substitution of the value of 
e from Eq. (4) gives 


€ 


4R/xc 


= ——_——._ rVcasin# 
1 + asi 
The outward normal force per unit arc is pV’. This is 
dN 4R/xc 


pV°rca sin 


1 + (4R/xc) 


Fic. 1. Nomenclature for ring airfoil. 
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The liftwise component of this normal force, expressed 
per unit span, is the same. Adding contributions from 
top to bottom introduces a factor of two: 

dL 4R/xc 


dy + (4R/mc) 


pV°xca sin 


The lift is seen to be distributed elliptically, and by 
integration the total value is 


4R/ xe 
i (5) 
It will be convenient to introduce a reference elliptic 
flat plate of span Rand area S’ = (x/2)Rc so that the as- 
pect ratio is A’ = 8R/xrc. Substituting in Eq. (5) gives 


Cy’ = L/(*/2)pV?S’ = 4n[A’/(A’ + 2)]a (6) 


Eq. (6) gives precisely twice the lifting-line result for 
the lift of the flat reference wing. Thus, the lift of the 
ring airfoil is twice the lift of an elliptic flat plate that 
spans a diameter and has a quarter the area. Also, the 
angle of downwash and the ratio of lift to induced drag 
are the same for the two airfoils, as can be seen either 
from Eq. (4) or from considerations of momentum. 

Eq. (5) or (6) for the lift applies to ring airfoils with 
diameter large relative to the chord. It may be expected 
to break down at very low aspect ratios A’ = 8R/xc 
just as the corresponding equation for flat elliptic wings 
breaks down. At these very low aspect ratios the wake 
would align itself with the axis of the ring, just as it 
does with the plane of a flat airfoil. The lift of both 
would then be independent of aspect ratio and would be 
determined only by a and the span (or a and the radius) 
(see references 3 and 4). From considerations of 
momentum the lift of a very low aspect ratio ring 
should be 
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L = 2prR?V?a 


In Eq. (5) if ¢ is allowed to become very large relative 
to R, 


L = 4prR?V*a 


or twice the preceding value. The lift of a very low 
aspect ratio flat plate of semispan = R should be 


L = prR?V?a 


The lifting-line theory for the elliptic flat plate 
gives 


L = [("/2)(o \(24a)[A’/(A’ + 2)] 
and for very low aspect ratio (A’ — 0) this becomes 
L = 
= 


again twice the correct value. Thus, even at very low 
aspect ratios, when Eqs. (5) and (6) are invalid, the lift 
of the ring airfoil is still twice the lift of an elliptic flat 
plate that spans a diameter and has a quarter the 
area. 
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Face Buckling and Core Strength Require- 
ments in Sandwich Construction’ 


CONRAD C. WAN} 
Chance Vought Aircraft Division of United Aircraft Corporation 


SUMMARY 


Relations between the wrinkling stress of the facings of sand- 
wich materials and the physical properties of the core are 
studied. The effect of initial waviness on the wrinkling stress 
and the strength requirements of the core are treated by the in- 
troduction of a numerical coefficient K. Charts for estimating 
wrinkling stresses are presented. 


INTRODUCTION 


ee HAS BEEN considerable interest during the 
past few years in the possible use of sandwich ma- 
terials for the primary structures of a modern airplane. 
Such sandwich materials consist of two thin layers of 
high-strength facings bonded to a thick layer of low- 
density core. The core, which provides shear connec- 
tion between the facings, increases the effective moment 
of inertia of the facing cross section and thus makes 
possible a considerable improvement in the stability 
of structural members. Suitable core materials can 
be found to perform the additional fuyction of stabiliz- 
ing the facings so that they will not wrinkle until a 
high state of stress in the material is reached. It is 
apparent that the core plays an important role in sand- 
wich materials, and it is desirable to know the strength 
and stiffness requirements for an efficient core material. 

The object of this paper is to derive certain theoreti- 
cal relationships between the wrinkling stress of the 
facings and the physical properties of the core in sand- 
wich construction so that a proper core material may be 
selected by means of these criteria. Surface waviness 
is incorporated in the present analysis, since it has 
been found that this waviness may have a detrimental 
effect on the wrinkling stress. Sandwiches under uni- 
form simple compression are treated first to establish 
the fundamental relationships. These are later ex- 
tended to include the case of sandwiches under biaxial 
stresses. 


GENERAL DISCUSSION 


The wrinkling phenomenon of the facings of sand- 
wich materials has been treated by various authors 
as a stability problem of thin sheets on an elastic foun- 
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dation. They have all made the usual assumptions 
that every component of the sandwich is homogeneous 
and perfect in its geometric configuration and that the 
core material possesses a finite stiffness but infinite 
strength. Consequently, their studies concern only 
classical buckling and yield the upper limit of the 
wrinkling stress. A brief review of their conclusions is 
given in the following. 


Gough, Elam, and deBruyne! and Williams, Leggett, 
and Hopkins’ treated with different degrees of rigor 
the case of a practical sandwich—i.e., a sandwich with 
a core of finite thickness with facings on both sides and 
complete adhesion between the faces and the core. The 
wrinkling stress of the facings was given in terms of a 
critical wave length, the properties of the facing ma- 
terial, and the Young’s modulus of the core material. 
Their analyses were mainly concerned with sandwiches 
of metal facings and onazote or plywood core. Balsa 
core was discussed by the second group of authors, but 
it was approximated by “‘rather artificial materials .. . . 
that only resemble the true materials insofar as their 
isotropic qualities allow.” They concluded that, for 
current designs, sandwich panels had no advantage 
over the more conventional sheet-stringer construction. 
The statement is believed to be due primarily to the 
fact that the balsa core has been treated as an isotropic 
material with a density of 10 lbs. per cu.ft., a shear 
modulus of 20,000 Ibs. per sq.in., and Poisson’s ratio 
of 0.3. The Young’s modulus of this arbitrary mate- 
rial has been calculated as 52,000 Ibs. per sq.in. Actu- 
ally, for balsa of this density the with-grain modulus is 
350,000 Ibs. per sq.in., while the cross-grain modulus is 
20,000 Ibs. per sq.in. It will be shown later that the 
Young’s modulus of the core in the direction perpen- 
dicular to the faces is critical and that considerable im- 
provement in the wrinkling stress can be anticipated 
if the core is arranged with the grain normal to the 
faces, thus taking advantage of the higher modulus in 
the with-grain direction. 


Hoff and Mautner* have recently reported work done 
on sandwiches of Papreg faces and cellular cellulose 
acetate core. They concluded that the wrinkling stress 
was independent of the geometric configuration of the 
sandwich and dependent only upon the moduli of the 
materials for the facings and the core. The following 
formula was recommended for a conservative estimate 
of the wrinkling stress: 


531 
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Ca, 1/, EEG, 


where 
Ger, = wrinkling stress in the facing 
E, = Young’s modulus of the facing 
E. = Young’s modulus of the core in the direction 


perpendicular to the facing 
G, = shear modulus of the core 


If the proper values are substituted in the above for- 
mula, the wrinkling stress for sandwiches of aluminum- 
alloy faces and 10 Ibs. per cu.ft. end-grain balsa (E, ~ 
10,000,000 Ibs. per sq.in., E, ~ 350,000 Ibs. per sq.in., 
G, = 20,000 Ibs. per sq.in.) will be 212,000 Ibs. per sq. 
in. Even when a reasonable reduced modulus is in- 
troduced, the calculated wrinkling stress is still well 
above that observed in practice. 

It is apparent that further study of theories concern- 
ing the wrinkling stress of sandwich materials has to be 
made so that a better understanding of the strength, 
as well as the stiffness, requirements of the core can be 
acquired and a closer correlation be obtained between 
the predicted strengths and those actually observed in 
test. In order to realize fully the high wrinkling strength 
obtainable in sandwich materials, it is desirable to have 
a core with a high strength and stiffness perpendicular 
to the facings, since the lateral restraint on the facings 
is due primarily to the core properties in that direction. 
Inasmuch as with wood and many of the synthetic core 
materials properties can be improved in one direction 
at the expense of those in the other two, sandwiches 
with aelotropic cores are discussed in detail under the 
section of ‘‘Analysis and Results,’’ while the results for 
isotropic cores are summarized in the Appendix. 
Wrinkling of the facings in the plastic range or the ef- 
fect of variable modulus is considered so that a solution 
can be obtained for any practical face material includ- 
ing metal, Fiberglas, Papreg, etc., provided the stress- 
strain relations are known. The Karmdn reduced 
modulus is used in this connection. Although this re- 
duced modulus applies only to struts with no initial 
curvature, it is believed that it offers the best means to 
express the plastic behavior of the facing material. 

It is the contention of the author that the wrinkling 
stress of the facings is closely related to the initial 
waviness present in the facings. It is a well-known fact 
that curved struts possess no definite buckling load, 
and the ultimate load they can sustain depends some- 
what on the yield strength of the material. Sandwich 
materials with initial waviness in the facings exhibit 
similar phenomenon. Most types of facings, especially 
metal, are made from commercial sheet stock, and as 
such they are usually not truly flat. Neither is the core 


material perfectly uniform in thickness, nor the material 
used for bonding. These initial irregularities may be 
corrected in part during the fabrication process, for 
example, by using stiff platens in the forming press, 
but the initial stresses set up in the bonding materials 
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act in much the same way as initial waviness when the 
sandwiches are out of the press. It is apparent that 
initial waviness is present and might plausibly be ex- 
pected to play a part in the wrinkling of the facings. 

The physical picture of the initial waviness can be 
further clarified by the following considerations. Let 
the facings of a sandwich be loaded in compression. 
Any wave that might be present in the facings in the 
direction of the load will be amplified. This tends to 
set up tensile forces in the bond because the facings are 
restrained by the core from bowing out (or peeling off). 
This amplification will become greater as the load is in- 
creased, and there will soon be a time when the ultimate 
tensile strength of the bond (or the core) is exceeded. 
The facings will then separate from the core, and the 
failure of the complete sandwich occurs. Similar rea- 
soning leads to the possibility of wrinkling of the facings 
because of failure of the core in compression or in shear 
rather than in tension. The theoretical wrinkling 
stress given by previous authors!~* may not be realized 
in all cases. 

The amount of initial waviness present in the facings 
can vary over a wide range. However, it is not im- 
possible to derive a rational expression to represent its 
general characteristics. It seems logical to assume that 
the amplitude of the initial waviness is greater for 
longer wave lengths and for thinner sheet thickness. 
The initial wave length can be taken to be the same as 
the wave length of the wrinkle in the facings, since at 
wrinkling (or failure) this is amplified to the greatest 
extent. The following expression, which has been in- 
troduced by Donnell‘ in the study of buckling of thin 
cylinders under axial compression, will be used in the 
present study: 


Wo (1) 
where 


Wo = amplitude of the initial waviness, in. 


X = wave length of the initial waviness, in. 

t,; = thickness of the facings, in. 

K = numerical constant—a measure of the wavi- 
ness 


It can easily be seen that this expression satisfies the 
characteristics described above. It has the additional 
advantage that any local weakness of bond or core can 
be incorporated in the empirical constant K when the 
wrinkling stress is to be evaluated. 


ANALYSIS AND RESULTS 


Studies have been made on sandwiches of identical 
facings with aelotropic or isotropic cores. The results 
from rigorous treatment of boundary conditions are 
extremely cumbersome for practical application. An 
approximate solution has been found for the case of a 
sandwich with end-grain balsa core which makes pos- 
sible the construction of a simple chart showing the re- 
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x 
te to 
Fic. 1. Forces acting on a typical sandwich under compression. 


lations between the various parameters associated with 
the wrinkling stress of the facings of a sandwich. This 
chart (Fig. 1) can be used for any core material exhibit- 
ing properties similar to end-grain balsa. This ap- 
proximate solution is presented in this section to illus- 
trate the effect of the initial waviness, while the more 
rigorous treatments are summarized in the Appendix. 

It is assumed that the core material under considera- 
tion has a predominant Young's modulus in the direc- 
tion perpendicular to the facings. The shear modulus 
and the Young’s modulus in the transverse direction are 
considered to be zero. This is practically the case for 
end-grain balsa where the shear modulus and the 
Young’s modulus in the transverse direction are, respec- 
tively, about one-twentieth and one-sixteenth of the 
Young’s modulus in the longitudinal direction. The 
sandwich is assumed to wrinkle in the symmetrical at- 
titude with a tension failure in the core. This assump- 
tion is based on numerous tests on sandwiches with 
aluminum-alloy faces and end-grain balsa core. It has 
been confirmed by comparing with the rigorous solution 
that, for balsa core, tension failure in the core occurs at 
practically the same facing stress for both the symmet- 
rical and the antisymmetrical attitude, while shear 
failure in the core corresponds to a much higher facing 
stress. 

The facings have the initial shape of 


wo = Wo sin ax (2) 
where 
Wo = [Eq. (1)] 
a = 2nr/r 
X = wave length of wrinkle 


When the facings are loaded in compression, the ad- 
ditional deflection shape can be taken as 


w= Wsin ax (3) 
The strain energy stored in a sandwich strip of unit 


width and one wave length is 


*/d>w\? 


2t,J 0 (uw) 2°L12(1 — te 


t, = thickness of facing, in. 


t, = thickness of core, in. 


TE, = reduced modulus of facing material, Ibs. per 
sq.in. 

uy = Poisson’s ratio of facing material 

E, = modulus of core-bond combination, Ibs. per 


sq-in. 


The virtual work done by the external forces (ot,) due 
to a virtual displacement dw = 61 sin ax is 


2 
= (2) f, + w+ iu) | 


| +m) dx = dota(W + 


At equilibrium the virtual work is equal to the change 
of the strain energy caused by the virtual displacements, 
and it shows that 
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Fic. 2. Wrinkling stress o.r, for sandwiches of end-grain balsa 
and 24S-T81 Alclad or 75ST Alclad facings. 
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Fic. 3. Wrinkling stress ocr, for sandwiches of end-grain balsa and facings with constant modulus. 


W = Qd./2E, (5) 


The facing stress is given by the following expression 
after Eqs. (1) and (5) are substituted into Eq. (4): 


all + = (72) + 


2E te 
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The maximum amount of deformation consistent with — \te (6) 


the ultimate tensile strength Q, of the core-bond com- 
bination is 
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The practical wrinkling stress is the minimum value 
of « as given by Eq. (6). The minimization is carried 
out with respect to at,, and the wrinkling stress is given 


by 
Ger, = V 1 — &) (7) 


where 


= V2rE,Ed)/3(1 — uf te 
= Ko /2Q, 


‘The wave length of the wrinkle at failure is accordingly 


Q 
| 


wee 
| 


Ner, = — (8) 


It is to be noted that the wrinkling stress decreases 
with an increasing value of K. When K is set equal to 
zero, the wrinkling stress given by Eq. (8) is comparable 
to the theoretical values found by previous authors.!' * 

Eq. (7) can be given in the form of charts such as 
shown in Figs. 2 and 3. It is possible to construct 
charts for any material by noting that the intercepts 
are 1/o.,, on the K/Q, axis and 2rE,/3(1 — uf)o.,.2 on 
the t,/E.ty axis. These results are readily obtained 
from rearranging Eq. (7) in polar form. 

The charts given in Figs. 2 and 3 for 24S-T81 Alclad 
and 75ST Alclad faces are good only for tension or com- 
pression failure in the core or the bond. However, it 
will be shown in the Appendix that shear failures in the 
core are not probable if the ultimate shear strength of 
the core is greater than 12 per cent of its ultimate ten- 
sile or compressive strength. This implies that, when 
new core materials are under consideration, the ulti- 
mate tensile or compressive strength perpendicular to 
the faces should be considered to be of greater signifi- 
cance than the ultimate shear strength in developing 
high wrinkling stresses in the faces. It is realized, of 
course, that in any practical sandwich panel subjected 
to transverse loading, a high shear strength in the core 
is highly desirable. 

It can be readily seen from the charts (Figs. 2, 3, and 
4) that the faces of a sandwich can be stabilized to a 
higher stress by using a core of higher stiffness. It also 
indicates that for the same amount of initial imperfec- 
tion the wrinkling stress can be increased by using core 
material of higher ultimate strength. 


APPLICATION OF RESULTS—VALUE OF K 


The results of the present study are given in the form 
of charts in Fig. 2 for facings of 24S-T81 Alclad and 
75ST Alclad sheets and core of end-grain balsa. Re- 
sults for facings with constant modulus and end-grain 
balsa are given in Fig. 3. The case of sandwiches with 
isotropic core and facings of constant modulus and 
t./ty of 40 is shown in Fig. 4. The application of these 
charts is best illustrated by an example. 

Suppose that a sandwich specimen is made of two 
0.012 gage 24S-T81 Alclad sheets and a 0.50-in. core 
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Wrinkling stress o,r, for sandwiches of isotropic core and 
facings with constant modulus when ¢./4 = 40. 


Fic. 4. 


of end-grain balsa, whose ultimate strength in tension 
is 800 Ibs. per sq.in. and whose Young’s modulus is 
330,000 Ibs. per sq.in. It follows from Fig. 2 that, for 
K = 0.010, the wrinkling stress would be 62,500 Ibs. per 
sq.in. Conversely, if the facings are to sustain a com- 
pressive stress of 62,500 Ibs. per sq.in., the ultimate 
tensile and/or compressive strength of the core ma- 
terial should be at least 800 Ibs. per sq.in. for K = 
0.010. 

The amplitude and the wave length of the initial 
waviness for the above example can be readily calcu- 
lated by means of Eqs. (8) and (1). 
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Fic. 5. Comparison of typical test results for sandwiches of end- 
grain balsa core and 75ST Alclad facings with analysis. 
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he. = 7 X 0.012 X 
V 0.76 X 0.284 X 10.5 X 10°/62.5 X 10? = 
0.23 in. 
Wo = 0.010 X (0.23)2/4x? X 0.012 = 0.0011 in. 


These values are considered to be well within expecta- 
tions. 

The value of K to be used in conjunction with the 
charts should be determined from statistical survey of 
tests on representative specimens. Any possible effect 
due to variations of bond strength and core strength 
can then be incorporated in the empirical constants. 
Typical test results for sandwiches of end-grain balsa 
and 75ST Alclad facings are shown in Fig. 5. Wrin- 
kling stresses according to Fig. 2 are also shown for K = 
0.010, 0.015, and 0.020. It appears that reasonable 
correlation between the tests and the present analysis 
can be obtained if K is taken to be 0.015. 


EXTENSION OF RESULTS TO SANDWICHES UNDER 
BIAXIAL STRESSES 


The approximate solution given previously for sand- 
wiches under simple uniform compression can be ex- 
tended to sandwiches under biaxial stresses. Let x and 


-y denote the directions of the principal stresses and 


assume that o; is greater than o,. The wave length of 
the wrinkles in the x- and y-directions will be denoted 
by A, and A,. The expression for Wo is assumed to be 
given by the following expression: 


Wo = (9) 


At equilibrium, the expression for the principal stress 


(“) 9 2E de ~2 
_ 1+ 12(1 — wf) \t + ty 
Q. ty m Yz 
(10) 
where 

n = 

= 

Ya * (1 + m*)/m?/dz 


The minimum value of ¢, corresponding to any stress 
ratio m is determined by setting the first derivative of 
expression (10) with respect to y, and m, respectively, 
equal to zero, These conditions, after rearrangements, 
are: 


= 120. — (*) + — &n) (11) 


by 
1 2 
n= (12) 


2 


2m? 


where 
= V2rEEd,/3(1 — uf)te 
1+ m’Koo_ 
m 20. 


The critical wrinkling stress and the wrinkling wave 
length are, respectively, 
1 + m? 
Cr, cr, = 1 + Em? — &m) (13) 
2\2 
«= arty (1 m?) TEy 
+ m?) 3(1 HP) 


Solution for sandwiches under biaxial stresses can 
be found only by successive approximation. The pres- 
ence of the reduced modulus has introduced a great 
deal of complication to the trial-and-error process. 
Since the above analysis is a crude approximate one, it 
seems rather unjustified to use its results for practical 
application. Therefore, the following alternative is 
suggested. 

For equal compression in both directions (n = 1), m 
must be equal to unity from consideration of sym- 
metry. This set of values satisfies Eq. (12). From 
Eqs. (13) and (14), the critical wrinkling stress and 
wave length are, accordingly, 


Ger, = Or = oy = V1 + 48% — 28) (18) 
Ne, = Az = Ay = — (16) 


It can be seen by comparing Eqs. (15) and (7) that the 
effect of equal biaxial compression is equivalent to 
doubling the value of the constant K. This indicates 
the possibility of applying a correction factor to K for 
biaxial stress conditions. The curve shown in Fig. 6 is 
believed to be on the conservative side for practical 
use. When the stress ratio is positive, the correction 
factor is assumed to vary linearly between the known 


1.0 -5 5 10 


Fie. 6, Correction factor to be applied to K for the determina- 
tion of ¢,, or, for biaxial loading. 
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end conditions. When the stress ratio is negative, the 
stiffening effect of the tensile stress is neglected. 


CONCLUSION 


An attempt has been made to study the strength re- 
quirements of core materials in sandwich construction. 
From consideration of the effect of initial waviness in 
the facings of sandwiches, a rational relation has been 
established between the ultimate core strength (or 
bond strength), the core modulus, and the wrinkling 
stress of the facings. A numerical constant K is in- 
troduced as a measure of the waviness of the facings. 
The magnitude of K for sandwiches with aluminum- 


alloy faces and end-grain balsa, as made with present. 


production technique, has been found to lie between 
0.01 and 0.02. Laboratory specimens might corre- 
spond to lower values of K. 

Charts have been presented for wrinkling stress of 
sandwiches with 24S-T81 Alclad and 75ST Alclad 
facings and end-grain balsa core. These charts can 
also be used for establishing strength requirements 
for core materials. 
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APPENDIX 


A rigorous solution of the problem considered in the 
present paper is given in the following. The funda- 
mental equations have been given by Williams, Leggett, 
and Hopkins.? The details of the derivation will be 
deleted for brevity. Necessary changes of notations 
have been made. 


ty = thickness of each facing 

i = thickness of core 

nN = wave length of wrinkle 

tE;, = Young’s modulus and Poisson's ratio of 
facing material 

th = Young’s modulus and Poisson's ratio of 
isotropic core 

E, = Young’s modulus of aelotropic core in 
the direction perpendicular to the fac- 
ings 

£,’ = Young's modulus of aelotropic core in 
the direction parallel to the facings 

G, = modulus of rigidity of aelotropic core 


v = E./G, 

Be = Poisson’s ratio of aelotropic core asso- 
ciated with EF, 

o = compressive stress in facing 

Sin. = wrinkling (or critical compressive) stress 
in facing 

Ue, We = displacements of points in the core 

Us, Wy, = displacements of points in the middle 


surface of the facing 

Or, O,, Tz; = components of stress in the core 

Q. = ultimate tensile or compressive strength 
of the core-bond combination, which- 
ever is lower 


S, = ultimate shear strength of the core 

a = 

B = at,/2 

6 = E,/E, 

2t,/ te 

K = measure of waviness in the facing 

Wo = Ka*/4r*t, 

Wo = initial shape of assumed waviness in the 
facing = Wo sin ax 

q, $ = values of o, and r;, at the common sur- 


face of the facing and the core 

These equations are: 

TE gt? d‘wy, d*(wy Wo) 
12(1 — uf’) dx dx 2 dx 

S = 
uy + (t,/2) (dw,/dx) = (A-3) 


and the stress-strain relations between ¢,, ¢:, 7;, and 
Ue, We. 

For isotropic cores, the stress distribution in the core 
can be described by stress function of the following type: 


@ = [A,az sinh az + A: cosh az + B, az cosh az + 
B; sinh az} sin ax (A-5) 


For aelotropic cores, the stress function can be as- 
sumed as 


o= [41 cosh + A: cosh z+ 
sinh — z + Bs sinh — sin ax (A-6) 
Qa) a2 


where a’, a:* are the roots of the equation 


A? — (A-7) 


Symmetrical and antisymmetrical attitudes of wrin- 
kling can be treated, respectively, by setting the B,'s 
and the A,’s equal to zero. 

All the pertinent expressions for uy, wy, g, s can be ob- 
tained by straightforward manipulations. Eqs. (A-1) 
and (A-2) are then reduced to the following system of 
equations: 


° 

K = 
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The values of A,’s and B,’s are determined, and they 
are introduced into the expressions for o, and 7;:. 
The magnitudes of these stresses are easily calculated 
at the face and the centerline of the core. Equating 
these stresses to the corresponding ultimate strength for shear failure. 
of the core, the stresses in the facings at equilibrium are The functions F(8, 7), G(8, r), and S(@, r) are sum- 
expressed by the following formulas: marized as follows: 


0 for tension failure. 


o KE, 
E. = F(g, T) + “5, (A-10) 


(A) Antisymmetrical Attitude 


(1) Isotropic Core: 


| 4— 
F(8, 7) = —8 + tanh + tanh? B + + (1 + (1 — tanh + 


61 
(1 + + be)? — we) (1 + Me) 
— (1 + tanh? ab + tanh? 3) +5 tanh 3] 
G(B, 7) = 28¢ + + — tanh? 8) + (3 — + tanh 8] 
Q(8, 7) = —8 + tanh + 6 tanh? 6 + 75 gape tanh 


Es He) (gags i 


tanh 6 + 8 — (1 + + la ah at) 


(2) Aelotropic Core: 


1 1 = 2 2 
F(8, 7) = —— tanh + — tanh 8 + 76 | = + (ax? He) tanh 
+ tanh + (ao? — a,*)8¢ tanh tanh + | + 64¢4 tanh 
ay ay 12(1 — ay 
(a? + te)? 
h 
1 — 
7) = = Bt + 76 tanh — B _ Car? + 8¢? tanh 
2 a2 ay 
1 1 ( c 
Q(8, 7) = —— tanh 8 + — tanh + 76 = be) tanh tanh 
a a2 2a, Qa) 
(a2? — a”) tanh tanh 
a) 
1 
S(8, r) = - (sech sech + 76 + yu.) sech + u,) sech Be? + 


(oe + u-) sech 2 tanh B — (a,? + u,) sech tanh 
2 (2) 


(B) Symmetrical Attitude The critical compressive stress in the facing naturally 
The functions F(8, 7), G(8, r), and Q(8, 7) can be has to be the minimum value of the Eqs. (A-9) and 
obtained by replacing “tanh” with “‘coth.”” The func- (A-10). This can be determined graphically in the 


tion S(8, 7) is not needed for the symmetrical attitude. following manner. Eqs. (A-9) and (A-10) are rear 


KEe og (A-9) 
+ 
Tt 
fe 
gra 
out 
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Q(8, 
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| 
| 
vie 


(A-9) 


(A-10) 


> sum- 


naturally 
A-9) and 
ly in the 


are reat 


ranged as 
, 
“™ 
1_ , K S(8,7) 


Various values of 6 are then substituted into the 
above equations, and the resulting relations are plotted 
onao-7 diagram. The intersection with the o—7 curve 
for the facing material locates the solutiom for o corre- 
sponding to a particular value of 8 and K/Q, or K/S,. 
The minimum value of o is obtained from a plot of o 
versus 8 for various values of K/Q, or K/S,. 

A comparison between the approximate solution 
given in the test and the rigorous solution presented 
above is shown in Fig. 7. It can be said that, for end- 
grain balsa, the tension failure will always happen first 
and that the approximate solution can be applied with- 
out any serious error. 

This can also be deduced by comparing the function 
Q(8, r) and S(8, 7). The minimum wrinkling stress 
corresponds to values of 8 greater than 3. For 8 = 
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Fic. 7. Comparison of results for a sandwich of 0.012 in. 
24S-T81 Alclad facings and 0.48 in. end-grain balsa (E. = 250,000 
Ibs. per sq.in., we = 0.50, x = 16.4, » = 22.6) from approximate 
and rigorous formulas. 


coth 8 is practically the same as tanh 8, and sech 8 is 
of the order of '/,0, and the results shown in Fig. 7 are 
to be expected. 
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Momentum Relations in Propulsive Ducts’ 


PHILIP RUDNICKt 
Vanderbilt Unwersity and The Johns Hopkins University 


SUMMARY 


This paper discusses the quantity 


A(p + pV?) 
considered as a function of position along a steady ‘‘unidimen- 
sional’’ compressible flow through a duct. (A = area of cross 


section; p = pressure; p = density; V = velocity.) The func- 
tion is called ‘‘stream thrust’’ and has the fundamental property 
that the difference between its values at any two sections of a 
stream gives the resultant of all normal and tangential forces act- 
ing on the duct walls between those sections. A description of the 
function is followed by some illustrations of its usefulness in con- 
nection with supersonic ram-jet design. 


NOTATION 

p = pressure, lb. (mass) per sq.ft. absolute 

p = density, slugs per cu.ft. 

k = gas constant, ft.-Ib. (force) per slug-°R. 

= absolute temperature, °R. 

Y = ratio of specific heat at constant pressure to that at 
constant volume 

A = cross-sectional area, sq.ft. 

V = flow velocity, ft. per sec. 

¢ = sonic velocity, ft. per sec. 

F = stream thrust, defined by Eq. (3), lb. (force) 

Sa = air-specific impulse, defined by Eq. (11), lb. (force)- 
sec. per Ib. (mass) 

Sr = fuel-specific impulse, defined by Eq. (11), Ib. (force)- 
sec. per Ib. (mass) 

F = mechanical force, lb. (force) 

g = acceleration of gravity, ft. per sec.? 

W = total mass flow rate, lb. (mass) per sec. 

a = rate of mass flow of air, Ib. (mass) per sec. 

4 = rate of mass flow of fuel, Ib. (mass) per sec. 

M = Mach Number (= V/c) 

# = defined by Eq. (15) 


g(M) = defined by Eq. (6) 
¥(x) = defined by Eq. (16) 
Numerical subscripts refer to specific cross sections of a 


duct. 
Literal subscripts have the following meanings: 


stagnation 
c = critical—i.e., when VJ = 1 
bd burner drag 
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INTRODUCTION 


NALYSIS OF THE OPERATION of a jet motor requires 
a new emphasis on momentum relations,! in con- 
trast to the usual treatment of heat engines in which 
energy considerations predominate. This paper de- 
scribes a working concept involving momentum which 
is applicable in any case of jet propulsion and which has 
proved useful in the practical design of ram-jet motors. 
The next two sections formally define and describe 
this quantity, here called ‘‘stream thrust.’ The final 
section undertakes to illustrate its range of usefulness, 
The illustrations are drawn from a context of super- 
sonic ram-jet analysis and design but do not in any 
way constitute a systematic treatment of this subject. 
They are intended only to exemplify an analytical ap- 
proach, which is considered to have some novelty if 
this may ever be said of purely Newtonian mechanics. 
This point of view has been developed and used by the 
propulsion group at the Applied Physics Laboratory. 
For contributions to the material in its present form, 
particular acknowledgments are due to J. E. Cook and 
T. Davis. 


DEFINITION OF STREAM THRUST 


The process to be considered in this section is the 
steady-state flow of a gas through a duct, taking full 
account of its compressibility. For simplicity the fol- 
lowing will be assumed: 

(a) The gas obeys the usual equation of state 


p/p = kT (1) 


and has a constant specific heat ratio y. 

(b) The duct wall is a surface of revolution. 

(c) Velocity and other flow variables are uniform 
over each cross section of the duct, and transverse com- 
ponents of velocity are negligible. Skin friction and 
shocks may occur, however. 

Under these circumstances, the resultant of all nor- 
mal and tangential forces exerted by the flowing gas on 
the walls of the duct included between any two cross 
sections, (1) and (2) in Fig. 1, is axial in direction and 
has the magnitude F2 given by 


Fo = Ao(pe + p2V2?) — Ai(pi + 
or more compactly 
Fy = — P, 2) 


where 
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MOMENTUM RELATIONS 


(2) 

Far 
4 

(a) 


(b) 


Fic. 1. Illustration of analogy between stream thrust in fluid 
flow and compressional forces in a rigid rod. 

(a) Fluid flow through duct. P; and P» are entrance and 
exit stream thrusts. f is total force exerted by the duct walls 
on the stream. Fy is the force on the duct walls equal and op- 
posite tof. See Eq. (2). 

(b) Rigid rod. F, is force exerted on rod at one end. ff is 
additional distributed force on rod. F» is force exerted by rod at 
the otherend. F\, F:, and f are analogous to P;, P2, and f in (a). 


P = A(p + pV?) (3) 


is regarded as a function of position along the stream, 
with values P; and P: at the corresponding sections (1) 
and (2). Fx is to be considered positive in Eq. (2) 
when directed from section (2) toward (1), regardless of 
the direction of stream flow. 

Eq. (3) shows the function P to be made up of two 
parts: Ap equals the pressure force transmitted across 
a cross section, and ApV? equals the time rate of mo- 
mentum flow across the same section. Eq. (2) is 
simply an expression of the fundamental mechanical 
principle that the time rate of creation of momentum 
within a region is equal to the resultant of all external 
forces acting across its boundaries. Considered in an- 
other way, Eq. (2) shows the force exerted by the stream 
on a duct wall to be measured directly by a change in 
the function P between the entrance and exit cross sec- 
tions. This relation is sufficiently useful to warrant 
giving the quantity P a name. It will be called the 
“thrust of the stream” or ‘‘stream thrust,’ by its anal- 
ogy to a compressive force or thrust transmitted along 
a rigid rod and caused to vary along the length of the 
tod by the existence of distributed external forces (see 
Fig. lb). It should be apparent that the restrictions 
assumed above on the character of the flow are not es- 
sential and that a much more general definition of stream 
thrust is possible. Such a definition is stated in the 
Appendix. 

The relation stated by Eq. (2) is generated by an 
impulse-momentum conservation law, which requires 
that the stream thrust be constant along the flow ex- 
cept as it is altered by external forces. Close parallels 


LSIVE DUCTS 


IN PROPU 


to this should be recognized in each of the following, pos- 
sibly more familiar, relations. 

(a) Conservation of mass leads to a mass-flow func- 
tion, pVA, which is constant along the stream except 
where matter is added or removed, as, for example, 
where fuel is injected into an air stream. 

(b) Conservation of energy requires a constant rate 
of energy transport along the stream, symbolized by 
constant stagnation temperature 7, except where this 
is altered by some energy source or sink, such as heat 
transfer across duct walls or combustion within a flow- 
ing gas. 

(c) Along a stream with 7, constant, the second 
law of thermodynamics requires entropy also to be con- 
stant, symbolized by constant stagnation pressure /,, 
except where this is altered by an irreversible process 
such as friction or shock. 

Thus, stream thrust may be classed with mass flow, 
stagnation temperature, and stagnation pressure; these 
four parameters derive their importance from their con- 
nection with the four conservation laws cited. 


DEPENDENCE OF STREAM THRUST ON IV, 7,, AND M/ 


Eq. (3) can easily be transformed to give the stream 
thrust P in terms of mass flow (IV), stagnation tem- 
perature (7°), and Mach Number (\/). The result 
may be written 


P = (4) 
where 
W +1)... +1)We 
4 vg 
| 
3 
(m) 
| 
2 
| 
| 
| 
M 
Fic. 2. Relative values of stream thrust as a function of 


Mach Number; for y = 1.2 and 1.4. See Eqs. (4), (5), and (6) 
for definitions. 4/7; and JW, satisfy Eq. (7) and are a possible 
pair of Mach Numbers before and after normal shock in air. 
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is plotted in Fig. 2. It has a minimum at J/ = 
1 and is normalized to unity at that point. This curve 
shows the relative variation of P with J/ along a stream 
in the common case in which W and 7, (and therefore 
P,.) are constant. Thus, one may consider the abscissa 
in Fig. 2 correlated with position along a stream flowing 
through a converging-diverging nozzle, from subsonic 
to supersonic velocity (or the reverse). The behavior 
of y, or P, which first decreases and then increases, re- 
flects the fact that the axial forces against the walls of 
the nozzle are in opposite directions on opposite sides 
of the throat. Some alternative forms of the function 
yg are shown in the Appendix. 

P. is the minimum value of the stream thrust, 
achieved when JJ = 1. Of the two expressions given 
for it in Eq. (5), the latter, involving the critical sonic 
speed c,, is somewhat more basic, since it remains valid 
in the case where y depends on temperature, provided 
y is given the value appropriate to the critical state. 
Eq. (6) and the other part of Eq. (5) are not so readily 
adapted to use for variable y. 

Eq. (4) implies some interesting facts that may not 
otherwise be immediately apparent. For a stream 
along which P, is constant (i.e., with mass flow and 
stagnation temperature constant), the force on the duct 
walls between any two sections depends only on the 
values of Mach Number at the two sections, regardless 
of the extent or character of shocks or other irreversibili- 
ties that may lie between. In particular, if entrance 
and exit Mach Numbers are equal, the net force on the 
intervening duct walls is zero, again regardless of ir- 
reversibilities. 


ILLUSTRATIVE APPLICATIONS OF THE STREAM THRUST 
FUNCTION 


(a) Normal Shock 


The relation between Mach Numbers ahead of and 
behind a normal shock is readily obtained by means of 
the function ¢ in Eq. (6) and Fig. 2. Stream thrust, 
W, and 7, must, of course, all be continuous across a 
normal shock. Hence,e 


will give the relation of the Mach Numbers J/; and M2 
before and after shock. Eq. (7) might be much simpli- 
fied algebraically, but, if used in connection with the 
graph in Fig. 2, it affords a simple way of relating Mach 
Numbers across a shock, as illustrated in the figure. 

In a similar way, the pressure ratio across a normal 
shock may be derived from a curve of y against x [see 
Eqs. (15) and (16) in the Appendix]. x; ‘and x2. must 
be related by 
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¥(x1) = (8) 
and 

po/ pr = X1/Xe (9) 
(since A; = A2), where subscripts 1 and 2 again refer 


to points just ahead of and behind a normal shock. 


(b) Thrust Against the Expanding Section af an Exit 

Nozzle 

The thrust* of a jet motor may be calculated from a 
difference in stream thrust for exhaust and intake, as in 
Eq. (2). Taking the intake conditions to be constant, 
an additive change in the stream thrust at the exhaust 
will produce the same additive change in thrust of the 
motor. 

If, for example, a nozzle with no expanding section 
(giving J = 1 at the exit) is to be compared with a 
deLaval nozzle with sufficient expansion (about 1.2:1) 
to give M = 1.5 at the exit, Fig. 2 [or Eq. (6)] shows 
that the gain in exit stream thrust is 6 per cent. The 
thrust of the motor [F2, in Eq. (2)] is increased by the 
same amount of added force but in a greater proportion, 
depending on the magnitude of the intake stream thrust 
[P; in Eq. (2)]. To complete a comparison of the two 
nozzles, the alteration in external or drag forces would, 
of course, also have to be taken into account. 


(c) Criterion of Quality of Combustor Performance 


It is sometimes desirable to state the thrust producing 
quality of an exhaust gas stream by means of an inten- 
sive property that is independent of the mass flow and 
also of the geometry of the exit nozzle. The stagnation 
temperature is a first approach to such a property. 
However, more satisfactory criteria, though not ideal, 
may be derived from the quantity P, of Eq. (5). If 
P,, is understood to refer to the stream formed by com- 
bustion of fuel at rate f in air at rate a—i.e., with 


W=a+f (10) 
then the quantities S, and S, defined by 
aS,=fS,;=P, (11) 


are called, respectively, air-specific impulse and fuel- 
specific impulse. They are functions of the tempera- 
ture and composition of the exhaust, hence depend on 
the ratio f/a and the completeness of combustion. 5, is 
of major significance when the design objective is 
maximum thrust from a limited air supply, as in the 
case of short-range flight of a ram-jet motor. 5S; is of 
importance as a criterion of fuel economy, as in a long- 
range ram-jet. It is analogous to the term “‘specific 
impulse’ as applied to rocket fuels. 


* “Thrust” is used here and in section (d) as the total force on 
internal duct surfaces, based on absolute pressures. For practi- 
cal purposes, only gage pressure should be used, so that thrust 
may vanish when the duct contains still air at one atmosphere. 
Shift to this latter basis can be made straightforwardly but, for 
simplicity, is omitted from the above discussions. 
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MOMENTUM ‘RELATIONS IN PROPULSIVE DUCTS 


Since both S, and S;are derived directly from P,, they 
are in the first instance performance criteria that 
characterize an exhaust stream with 1 = 1, separate 
knowledge of y and c, in Eq. (5) not being necessary. 
When an expanding nozzle is used to give M > 1 at the 
exit, the added thrust derived from the expansion de- 
pends on y as wellas P,. However, as already pointed 
out, this added thrust is comparatively small for 
moderate expansions; therefore, in such cases, S, and 
S, are still important criteria of thrust producing 
quality, separate knowledge of y being a secondary 
matter. 

As well as being theoretically appropriate, S, and Sy, 
are rather directly determinable by expeiiments of two 
types: 

(1) Direct thrust measurements made statically on 
a combustion chamber with exhaust at “= 1—i.e., 
with nonexpanding exit—operated with known rates 
of air and fuel flow. If air and fuel are introduced 
transversely, P; in Eq. (2) becomes zero, and the meas- 
ured thrust directly determines the exit stream thrust 

(2) Measurement of pressure at the inlet end of a 
cylindrical combustion chamber with the exhaust end 
unconstricted, asin Fig.3. One must also know air and 
fuel flow rates and the total aerodynamic drag force 
within the chamber. The relations in this case are 
indicated more fully in the next section, Eq. (12). 


(dq) Some Elementary Aspects of Ram-Jet Design 


A simple ram-jet motor is indicated schematically in 
Fig. 3. Section 2 divides the structure into its two 
principal parts: the forward portion, the diffuser, 
which receives incoming air during flight, reduces its 
velocity relative to the duct, and thereby develops an 
increased pressure; and the after section, the combus- 
tor, in which fuel is injected and burned to produce 
thrust. In the case here considered, the pressure p» at 
the diffuser exit is supposed sufficiently high to produce, 
at the combustor exit, sonic velocity (17; = 1) at a pres- 
sure above atmospheric (though necessarily much below 
p2). 

Suppose that specific values have been chosen for the 
flight speed and the maximum cross section A» (so that 
the aerodynamic drag of the motor may also be consid- 
ered roughly fixed). One may then proceed to some 
rough design considerations. 

The propulsive thrust* of the motor is, by Eq. (2), 
a difference between stream thrust at the exhaust and 
at the inlet. The latter will be determined by the rate 
of air intake, which for maximum thrust should, of 
course, be as small as allowed by other considerations. 
The exhaust stream thrust may be expressed as the 
stream thrust at section 2, diminished by the drag force 
existing beyond that point, 


* See footnote, page 542. 


: DIFFUSER | COMBUSTOR 
ap) (2) ¢3) 


Fic. 3. Schematic diagram of a simple ram-jet motor. 


Ps = Ax(p2 + p2V2?) — Fog (12) 


Typically, the first of the three terms in the right mem- 
ber, A2p2, is larger by an order of magnitude than either 
of the following two and thus roughly determines the 
available stream thrust P;. We may then say that 
p:—and, hence, P; approximately—is fixed by the flight 
velocity and diffuser efficiency. The P3 so fixed is an 
upper limit to the possible thrust derivable from this 
motor. 


On the other hand, P; must be related to the combus- 
tor performance by 


P; =aS, (13) 


This relation permits the choice of an air rate a consist- 
ent with available combustion quality S, and the de- 
sired streain thrust P; allowed by the diffuser. There 
is an upper limit to air rates at which the pressure recov- 
ery in the diffuser is reasonably complete; Eq. (13) 
fixes a corresponding lower limit for a suitable S,. If 
available values of S, exceed this limit, the motor 
geometry indicated in Fig. 3 is feasible. If maximum 
thrust is required, the largest S, and smallest a will be 
chosen; small S, and larger a may give better fuel econ- 
omy. The intake area A; is adjusted to be consistent 
with the desired air rate. Many other design questions 
can be treated by similar considerations. 


APPENDIX 
(a) General Definition of Stream Thrust 


A generalized definition of a stream thrust vector can 
easily be written with the aid of dyadic notation. Let 
S be the stress dyadic, so defined that S-dA is the force 
transmitted from the negative to the positive side of 
the elementary area dA in the fluid (by pressure and 
viscosity). Then the stream thrust vector P is defined 
as 


P= + pVV)-dA 


The integral is taken over an arbitrary entrance or exit 
area A, with dA in the downstream direction. Eq. (2) 
may then be used as a vector equation for a duct of any 
shape. 
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(b) Alternative Forms for ¢(M) The function g(/) can also be expressed in terms of 
The simplest mathematical form is assumed by the VY: 
defined by Eq. it is expressed in — Dx? 
erms of the ratio The result is = = — G?— (16) 
y? — 1)x 
¢g(M) = [1 + (V/c,)?]/2(V/e,) (14) 


A curve of ¥(x) plotted against would be somewhat 
It is sometimes useful to have the function g(J/) ex- similar to Fig. 2, since x increases monotonically with 
pressed in terms of still another independent variable Vand is equal to M when both are equal to 1. 
namely, the product of pressure and area of cross sec- 
tion. A dimensionless parameter x involving this prod- RErEeence 
uct is defined by 


1 See, for example, Moss, Sanford A., and Foote, W. R,, 


ge 2 1 1/2 Some Elementary Details of Ramming Intake and Propulsive Jets for 


ite (y+ 1)Ap ? ++ 1 Airplanes, Journal of the Aeronautical Sciences, Vol. 13, No. 3, 
pp. 111-118, March, 1946 (a discussion of momentum relations 


(15) _ for incompressible flow). 
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Errata 


The following corrections to the article “Stress Analysis of Rotor Blades” (JouRNAL, Vol. 14, No. 6, p. 315, June 
1947) were submitted by the author, Gabriel Horvay. 

(1) Page 315.—Egs. (39a) and (39b) replace ¢ by y. 

(2) Page 317.—Under heading ‘‘Outline of the Tabular Procedure..... ”’ change third line to: “‘....; inflow 
velocity and lift curve slope are constant through the rotor disc; flapping hinge is uncocked and coincides with the 
center of rotation, higher than second harmonic effects can be neglected,” etc. 

(3) Page 319.—In Eq. (14d), subscript & belongs outside of the brace. 

(4) Page 320.—Eg. (18) should read: 


te = dx + dx?}, 


(5) Page 321.—Edq. (24) should read: 


6J = 0 etc.) 


(6) Page 324.—Eq. (42) should read: 


g=T/x=B=.... 
(7) Page 326.—The last line in the right column should read: ‘and the numerical co‘umns cos 0, cos 2y, sin 2y of 
Table 4b....” 
(8) Page 330.—In Eq. (44*), integrate from zero. 
(9) Page 330.—In the right column, the seventh line from the top should read: “using the deflection curve y = 


1 — cos !/orx, and making.... 
(10) Page 331.—Replace Eq. (71) by: 


tip 


(11) Page 331.—Add the following at the end. of the sentence following Eq. (71): In performing the calculations, 
it is convenient to replace the steps S- and S, by 


Si = —{Ce}, (71a) 


and 
tip 
Se = — (71b) 
k 
Evidently, 
Si + Se = So + Sy (71’) 


(12) Page 331.—In Eq. (75), double dot ¢ in the first term. 

(13) “Page 331.—In Eq. (80), double dot y in the last term of the first member. 

(14) Page 331.—In Eq. (81b), precede sin 2y by factor T,. 

(15) Page 331.—In Eq. (82b), replace by (¢), and add footnote: e.g., = 2wAs. 
(16) Page 331.—In Eq. (83), add subscript » to T. 

(17) Page 332.—In Eq. (84), replace c(x) by C(x). 
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